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- — r-£>“  The  SRI  International  four-vsvelength  (0.53,  1.06,  3.8,  10.6  ym)  Ildar 
aystms  mas  used  during  the  SNOW-ONE -B  and  Saoke  Week  XI/ SNOW-TWO  field 
experiments  to  validate  Its  capabilities  in  assessing  obscurant  optical  and 
physical  properties.  The  Ildar  viewed  along  a  horizontal  patL  terminated  by 
a  passive  reflector.  Data  examples  vers  analysed  in  tsrms  of  tlme-dspsndsnt 
transmission,  wavelength  dependence  of  optical  dspth,  and  range-resolved 
extinction  coefficients.  Thrss  methods  wars  ussd  to  dsrlys  extinction  data'- 
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from  the  lidar  signatures.  These  were  the  target  method,  Klett  method  and 
experimental  data  method.  The  results  of  the  field  and  analysis  programs 
are  reported  in  the  journal  and  conference  papers  that  are  appended  to  this 
report,  and  Include:  % 


•  ^Comparison  Study  of  Lidar  Extinction  Methods,4*  submitted  to 

Applied  Optics  ^ 

•  "Error  Analysis  of  Lidar  Solution  Techniques  for  Range-Resolved 
Extinction  Coefficients  Based  on  Observational  Data,"  Smoke/ 

i  Obscurants  Symposium  IX  , 

•  "Four-Wavelength  Lidar  Measurements  from  Smoke  Week  VI/SNOW-TWO," 

\  Smoke/Obscurants  Symposium  VIII 

>  ' 

*'  ""SNOW-ONE- B  Multiple- Wavelength  Lidar  Measurements,*"  Snow  Symposium 

;  iii  . 


**  *Lidar  Applications  for  Obscurant  Evaluations,4*  Smoke/Obscurrmts 
Symposium  VII. 

The  report  also  provides  a  summary  of  background  work  leading  to  this  project 
and  of  project  results. 
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ABSTRACT 


The  SRI  International  four-wavelength  (0.53.  1.06,  3.8,  10.6  um) 
lidar  system  was  used  during  the  SN0W-0NE-B  and  Smoke  Week  XI/SN0W-TW0 
field  experiments  to  validate  its  capabilities  in  assessing  obscurant 
optical  and  physical  properties.  The  lidar  viewed  along  a  horizontal 
path  terminated  by  a  passive  reflector.  Data  examples  were  analyzed  in 
terms  of  time-dependent  transmission,  wavelength  dependence  of  optical 
depth*  and  range-resolved  extinction  coefficients.  Three  methods  were 
used  to  derive  extinction  data  from  the  lidar  signatures.  These  were 
the  target  method,  Klett  method  and  experimental  data  method.  The 
results  of  the  field  and  analysis  programs  are  reported  in  the  journal 
and  conference  papers  that  are  appended  to  this  report,  and  include: 

•  "Comparison  Study  of  Lidar  Extinction  Methods,”  submitted  to 
Applied  Optica 

•  "Error  Analysis  of  Lidar  Solution  Techniques  for  Range-Resolved 
Extinction  Coefficients  Based  on  Observational  Data,” 
Smoke/Obscurants  Symposium  IX 

e  "Four-Wavelength  Lidar  Measurements  from  Smoke  Week  VI/SNOW- 
TWO,"  Smoke/Obscurants  Symposium  VIII 

e  "SNOW-OME-B  Multiple-Wavelength  Lidar  Measurements,”  Snow 
Symposium  III 

e  "Lidar  Applications  for  Obscurant  Evaluations Smoke/Obscurants 
Symposium  VII. 

.The  report  also  provides  a  summary  of  background  work  leading  to  this 
project,  and  of  project  results. 


I  BACKGROUND 


Earlier  studies  conducted  for  ARO  used  single  (0.69  pm)  and  dual 
(0.69  and  10.6  pm)  wavelength  lidar  systems  to  observe  obscurant 
aerosols.  Typically,  the  liaar  viewed  along  a  horizontal  path  termi¬ 
nated  by  a  passive  reflector  with  a  uniform  surface.  The:  lidar  target 
returns  were  interpreted  in  terms  of  transmission  of  laser  ene:  ^y  as  a 
function  of  event  time.  The  transmission  data  was  then  correlated  with 
the  backscattered  signal  received  from  the  obscurant  aerosols.  Experi¬ 
ments  with  the  two-wavelength  lidar  during  field  tests  DIRT-1  and  Smoke 
Week  II  showed  that  the  wavelength  dependence  of  transmission  (extinc¬ 
tion)  and  backscatter  could  be  used  to  Identify  the  obscurant  as  smoke 
o'*  dust,  based  on  a  strong  wavelength  dependence  on  particle  size. 
However,  relatively  large  data  scatter  occurred  in  plots  of  optical 
depth  evaluated  from  the  target  returns  plotted  against  path-integrated 
backscatter.  Assuming  the  data  scatter  was  a  result  of  changing 
particle  sizes,  theoretical  considerations  indicated  that  a  3“to-4  pm 
wavelength  lidar  may  provide  higher  correlation  than  a  0.69  or  10.6  pm 
wavelength  system.  In  addition,  a  shorter  wavelength  would  be  useful  to 
determine  the  mean  particle  size  of  smoke  aerosols. 

A  four  wavelength  lidar  system  was  constructed  and  initially  tested 
at  SRI  under  ARO  Contract  DAAG-20-80-C-0003.  The  lidar  was  designed  to 
provide  information  on  propagation  at  wavelengths  commonly  used  by  DoD 
laser  systems,  and  to  provide  for  multiple  wavelength  analysis  in  terms 
of  particle  characteristics  and  aerosol  densities.  Pulses  at  0.53», 

1.06,  3.8*  and  i 0.6  pm  can  be  transmitted  along  the  same  path  within  a 
150  ps  interval  to  overcome  space  and  time  separation  problems  associ¬ 
ated  with  the  earlier  dual  wavelength  system.  Initial  tests  at  SRI 
indicated  the  lidar  was  ready  for  operation  during  scheduled  Army  Field 
tests  which  could  provide  the  supplementary  data  needed  to  validate  the 
lidar  performance  and  data  analysis  procedures.  This  project  provided 
that  opportunity. 


II  SUMMARY  OF  PROJECT'  RESULTS 


The  results  and  findings  obtained  on  this  project  are  reported;  in 
the  Journal  and  conference  papers  reproduced  as  Appendixes  to  this 
report.  A  brief  summary  is  presented  below. 

The  four-wave length  lidar  system  was  transported  to  Grayling, 

Michigan,  to  participate  in  the  SNOW-ONE-B  field  test  conducted  during 
November  and  December  1983.  Prior  to  this  time,  a  new  digital  data 
acquisition  system  was  installed  within  the  lidar  to  improve  digitiza¬ 
tion  accuracy  and  reliability.  The  lidar  operated  on  all  SNGW-ONE-B 
tests  without  failure.  Because  of  eye-safety  considerations,'  only  the 
1.06,  3*8  and  10.6  pm  wavelength  energy  were  transmitted.  The  lidar  van 
was  positioned  alongside  transmissometer  instrumentation  vans  operated 
by  the  U.  S.  Navy  (NRL)  and  the  U.  S.  Army  Atmospheric  Science  Labora¬ 
tory.  The  lidar  and  transmisslometers  were  aligned  to  observe  parallel 
paths  extending  over  a  600  m  range.  A  passive  reflector  (target)  was 
placed  at  the  far  end  of  the  lidar  path  to  provide  for  a  constant  re¬ 
flectivity  surface.  The  primary  objective  of  the  SN0W-0NE-B  tests  were 
to  characterize  electromagnetic  propagation  through  dense  snow.  Only 
light  snowfall  occurred  during  the  observational  period. 

Lidar  target  returns  were  interpreted  in  terms  of  transmission,  and 
results  were  compared  to  data  provided  by  the  NRL  transmisslometers. 
Excellent  agreement  was  obtained  as  shown  by  the  data  example  presented 
in  Appendix  D  collected  during  an  IR  screener  smoke  test.  The  lidar 

results  showed  that  the  smoke  attenuation  (optical  depth)  was  within_2Q) _ 

of  bein?  equal  at  the  1.06,  3.8  and  10.6  pm  wavelengths.  Data  collected 
in  snowfall  showed  that  optical  depth  at  10.6  pm  was  76)  greater  than  at 
1.06  pm,  and  59)  greater  than  at  3.8  pm.  This  enhanced  attenuation  at 
longer  wavelengths  is  probaLly  a  result  of  multiple  scattering  effects 
from  large  scattering  particles.  Much  of  the  large  particle  scattering 
is  in  the  forward  direction  and,  therefore,  more  scattered  energy  is 
directed  back  into  the  initial  laser  beam  at  longer  wavelengths  and 
results  in  less  attenuation. 

Plots  of  optical  depth  (target  returns)  against  Inter -ated  back- 


scatter  (aerosol  returns)  show  more  data  scatter  than  was  anticipated — 
especially  at  3.8  ym  where  theory  had  indicated  that  a  better  correla¬ 
tion  might  be  obtained.  However,  further  inspection  of  the  lidar  back- 
scatter  signatures  indicated  that  the  dynamic  range  of  the  log  ampli¬ 
fiers  used  in  the  lidar  receiver  was  not  sufficiently  large  to  accommo¬ 
date  the  strong  signals  returned  from  dense  aerosols,  snow,  and  the 
target.  Subsequently  the  four  decade  amplifiers  were  replaced  with  six- 
decade  units.  The  SN0W-8-0NE  data  was  also  used,  to  conduct  a  'hole* 
and  'patch'  analysis  of  smoke  distributions.  These  results  are  further 
discussed  in  Appendices  C  and  D. 

As  a  second  phase  of  the  prelect,  the  improved  four-wavelength 
lidar  system  participated  in  the  Smoke  Week  XI/SNOW-TWO  field  tests 
conducted  at  Grayling,  Michigan,  during  the  1983-84  winter.  The  lidar 
procedures  were  nearly  identical  to  those  used  during  the  SNOW-ONE-B 
tests  with  the  exception  that  permission  was  granted  to  transmit  laser 
energy  at  all  four  wavelengths.  The  lidar  operated  without  failure  on 
25  obscurant  tests  scheduled  between  10  and  16  January.  Unfortunately, 
a  guy  wire  that  intersected  the  lidar/target  path  was  installed  Just 
prior  to  the  test  series,  and  reflections  from  the  wire  Interfered  with 
aerosol  scattering  on  some  tests.  However,  the  lidar  signatures  showed 
more  variability  with  smoke  density  and  greater  range  response  resulting 
from  the  enhanced  dynamic  range  receivers.  Plots  of  optical  depth 
against  path- integrated  backscatter  had  substantially  less  data  scatter, 
indicating  that  the  lidar  response  to  aerosol  density  had  been  improved. 

The  analysis  program  concentrated  on  methods  to  derive  range- 
re.  olved  extinction  coefficients  from  the  lidar  signatures.  Three 
methods  were  used.  The  target  method  uses  the  target-derived  trans¬ 
mission  and  path-integrated  backscatter  to  evaluate  a  single  Calibration 
constant  for  each  lidar  observation.  The  calibration  constant  can  then 
be  used  with  integrated  backscatter  to  evaluate  extinction  along  the 
lidar-to-target  path.  This  method  is  based  on  the  single-scatter  lidar 
equation,  and  provides  the  best  estimates  of  extinction  along  the  path 
when  multiple  scattering  is  not  Important.  The  method  does  provide  some 
correction  for  multiple  scattering,  as  the  constant  is  evaluated  from 
the  actual  lidar  signature.  It  was  shown  that  the  solution  constant 


varies  from  observation-to-observation  so  that  a  universal  constant 
cannot  be  u?ea  to  evaluate  single-ended  lidar  observations. 

The  second  method  (Klett  method)  ‘  ;  a  modification  of  a  formal 
solution  to  the  single-scattering  >  'dar  equation.  The  solution  assumes 
a  known  relationship  between  backscatter  and  extinction  coefficients  and 
extinction  near  the  far  end  of  the  lidar  observation.  Deviations  (RMS 
differences)  between  the  Klett  and  target  method  were  computed  for 
various  solution  parameters.  The  resulting  extinction  profile  for  solu¬ 
tion  parameters  that  yielded  minimum  RMS  differences  agreed  well  with 
the  profile  derived  from  the  target  method.  This  was  expected,  as  both, 
methods  are  derived  from  the  single-scattering  lidar  equation.  Contour 
plots  of  RMS  values  as  a  function  of  the  two  solution  parameters  were 
constructed,  as  were  contour  plots  of  transmission  differences.  These 
plots  show  that  the  Klett  method  is  relatively  sensitive  to  expected 
errors  in  assumed  solution  parameters. 

The  third  method  was  the  experimental  data  method  which  established 
an  experimental  relationship  between  the  optical  depth  and  the  path- 
integrated  backscatter  evaluated  from  the  lidar  signature.  Therefore, 
this  method  accounts  for  multiple  scattering  effects.  Comparison  of 
extinction  profiles  with  those  derived  from  the  target  method  showed 
that  the  experimental  data  method  provided  for  more  attenuation  In  the 
near  side  of  an  obscurant  cloud  and  for  less  attenuation  in  the  far 
side,  indicating  the  importance  of  multiple  scattering.  Before  any  one 
method  is  Judged  superior,  an  experiment  shpuld  be  conducted  to  validate 
the  lidar  methods  for  both  low-  and  high-density  obscurant  events. 

Finally,  the  range-resolved  extinction  coefficients  evaluated  as  a 
function  of  obscurant  event  time  were  used  to  construct  a  contour  plot 
of  the  time  and  range  variability  of  smoke  distribution  along  the  pro¬ 
pagation  path.  Such  data  are  useful  in  evaluating  obscurant  effects  on 
electrooptical  systems  operating  within  the  obscurant  cloud. 

A  more  detailed  accounting  of  project  results  and  data  examples  is 
presented  in  the  Appendices. 


4 


•jM-, 


Ill  RECOMMENDATIONS 


The  liaar  technique  is  the  test  method  available  for  evaluating 
spatial  and1  temporal  distributions  of  obscurant  properties  and  obscurant 
effects  on  military  laser-based  weapon  systems.  However,  uncertainties 
in  relating  backscatter  data  to  extinction  and  density  remain. 
Multiple-wavelength  systems  provide  additional  information  to  help 
analyze  lidar  signatures  in  terms  of  needed  obscurant  parameters. 
However,  this  approach  is  limited  because  of  different  1 idar 
characteristics  of  multiple  laser  and  detector  systems,  and  because 
newer  military  obscurants  do  not  follow  available  scattering  theories. 
Nevertheless,  the  usefulness  of  lidar  for  obscurant  evaluations  nas  been 
demonstrated  and  further  development  is  warranted. 

The  following  recommendations  are  offered  to  further  develop  the 
lidar  technique  and  enhance  its  value  to  military  programs. 


•  Coordinated  in  situ  sensing  of  aerosol  densities  along  lidar 
propagation  paths  to  provide  data  needed  for  validation  and/or 
development  of  lidar  signature  analysis  methods. 

•  Development  of  methods  to  deconvolute  lidar  performance  factors 
from  lidar  signatures  so  that  multiple  wavelength  data  can  be 
analyzed  effectively  in  terms  of  atmospheric  optical  and 
physical  properties. 

•  Addition  of  vertical  scanning  optics  to  the  four-wavelength 
lidar  system  to  investigate  obscurant  distributions  in  the 
vertical  and  their  dependence  on  meteorological  conditions  and 
terrain  properties. 

•  Investigation  of  airborne  lidar  systems  for  measurement  of  wide- 
area  screening  smeke  distributions, 

•  Continued  participation  in  the  Smoke  Week  series  and  reporting 
of,  results  at  Smoke/Obscurants  Symposia  to  promote 
interrelationships  between  lidar  development,  obscurant 
evaluations  and  electroopt,cal  weapons  testing  programs. 
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COMPARISON  STUDY  OF  LIDU?  EXTINCTION  METHODS 


E.  E.  Uthe  and  J.  M.  Livingston 
SRI  International 
Menlo  Park,  CA  9^025,  USA 


ABSTRACT 


The  SRI  four-wavelength  (0.53,  1.06,  3.8  and  10.6  yin)  Ildar  system 
was  used  during  the  Saoke  Week  XI /SNOW-TWO  field  experiments  to  observe 
range-resolved  backscatter  and  target-raflected  laser  energy  along  a 
horizontal  path  Intersected  by  test  obscurant  aerosols. 

Target-reflected  laser  energy  provided  an  excellent  measure  of 
time-dependent  path-integrated  extinction  (e.g.,  transmission) .  These 
values  allow  the  range-resolved  backscatter  to  be  analyzed  in  terms  of 
range-resolved  extinction  when  single-scattering  is  assumed.  This  two- 
ended  (lidar  and  target)  method  provides  extinction  coefficients  along 
the  optical  path  that  are  consistent  with  Independent  transmlssometer 
observations.  Tne  values  are  the  most  accv*ate  that  can  be  derived  from 
the  lidar  technique  when  multiple  scattering  is  limited. 


Extinction  profiles  evaluated  from  the  lidar/ target  method  can  be 
utilized  to  investigate  the  limitations  and  accuracies  of  algorithms 
which  are  used  to  derive  extinction  data  from  lidar  backscatter 
signatures  when  no  target  return  is.  available.  Single-ended  lidar 
methods  are  needed  for  vertical  and  slant-path  measurements.  Two 
single-ended  techniques  were  investigated.  The  first  is  the  well-known 
Klett  modification  of  the  analytic  solution  to  the  single-scattering 
lidar  equation.  The  Klett  technique  can  reproduce  the  extinction 
profile  derived  from  the  lidar/target  method  as  expected  since  both  are 
derived  from  the  single-scattering  lidar  equation.  However,  this 
requires  relatively  accurate  information  on  a  boundary  value  of 
extinction  and  the  backseat  ter-to-eati  net  ion  relationship.'  Contour 
plots  are  presented  which  show  deviation  of  the  Klett  extinction  profile 
from  the  lidar/target  extinction  profile  as  a  function  of  the  solution 


parameters. 

A  seoond  method  of  deriving  extinction  profiles  from  single-ended 
lidar  measirements  is  based  on  an  observed  relationship  of  optical  depth 
to  path- Integrated  lidar  algnal.  Such  a  relationship  derived  from 
lidar/target  data  Includes  multiple  scattering  effects  on  the  extinction 
evaluation.  The  method  does  not  reproduce  the  extinction  profile 
derived  from  the  aingle-eoattering  lidar/target  and  Klett  methods  — 
indicating  the  importance  of  multiple  scattering  on  the  lidar 
baoksoatter  signature. 
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I  LIDAR  DATA  COLLECTION 


The  SRI  four-wavelength  (0.53.  1.06,  3*8  and  10.6  urn)  Ildar  system 
was  used  during  the  Smoke  Week  XI/SI.'OW-TWO  field  experiments  to  collect 
data  for  analysis  of  obscurant  extinction  distributions.  Basically, 
laser  pulses  at  the  four  wavelengths  were  transmitted  along  a  horizontal 
path  terminated  at  a  distance  of  600  m  from  the  lidar  by  a  target  of 
uniform  reflective  properties. 

The  laser  pulses  at  each  wavelength  were  emitted  sequentially  along  the 
same  path  within  a  60  us  interval, to  reduce  effects  of  time-  and  space- 
varying  aerosol  concentrations  on  each  four-wavelength  pulse  set.  About 
60  pulse  sets  were  transmitted  for  each  minute  of  operation.  Both 
energy  backscattered  from  atmospheric  aerosols  and  energy  reflected  from 
the  solid  target  was  collected  within  a  telescope  and  directed  to 
appropriate  detectors  using  dichroic  filters.  After  logarithmic 
amplification,  backscatter  signatures  were  digitized  by  an  8-bit  100-MHz 
transient  recorder,  and  the  resulting  digital  records  were  stored  on 
nine-track  magnetic  tape  for  subsequent  analysis.  Obscurant  aerosols 
were  generated  upwind  of  the  lldar/target  path,  and  lidar  observations 
were  Initiated  Just  before  the  obscurants  .  Intersected  the  viewing 
path.  Other  details  of  the  lidar  system  and  its  application  to  the 
Smoke  Week  VI/SNOW-TWO  field  test  have  been  presented  elsewhere  (JJ . 

II  DATA  ANALYSIS  METHODS  AND  RESULTS 

A.  Target  Method 

Target  returns  observed  by  a  lidar  system  can  be  Interpreted  in 
terms  of  transmission  across  the  lldar-to-target  ^th.  Previous  results 
have  shown  that  lidar  derived  transmissions  agree  with  Independent  data 
records  collected  with  conventional  two-ended  transmissometers  (2). 
Therefore,  the  four-wavelength  lidar  target  returns  collected  C'—tng  an 
obaourant  event  provide  the  time  history  of  transmission  a';  ec^h 


wavelength,  and  provide  data  to  evaluate  optical  depth  (u--ln  T) 
relationships  between  wavelengths,  as  shown  by  the  example  presented  in 
figure  1 . 

Fernald  et  al.  (3),  among  others,  have  shown  that  if  the  ratio  of 
the  extinction  coefficient  to  the  backscatter  coefficient,  k«o(r )/3(r ) , 
remains  constant  over  the  range  interval  Ar  being  investigated,  then  the 
two-way  transmittance,  T^(r),  follows  directly  from  differentiation  of 
the  basic  lidar  equation,  and  is  given  by 

T2(R)  -  1  -  ~  J  P(r)dr  (1) 

0 

In  this  equation,  P(r)  is  the  signal  whose  amplitude  is  proportional  to 
the  power  received  from  a  scattering  volume  at  range  r,  and  is  corrected 
for  the  geometric  range-squared  effect  and  .  transmitted  energy 
variations.  C  is  a  Ildar  calibration  constant. 

The  constant  k/C  can  be  determined  by  combining  the  known  two-way 
transmittance  between  the  laser  and  the  target  with  the  Ildar  return 
signature  integrated  over  the  same. range.  The  value  of  the  constant  k/c 
may  then  be  substituted  back  into  (1),  and  combined  with  the  integrals 
of  the  received  backscatter  signals  to  Intermediate  ranges  to  yield  a 
horizontal  transmission •  profile.  Since  transmission  is  directly 
related  to  optical  depth  through  T2(R)  -  e-2u^,  where  '  optical 
depth  u(R)  •  J  o(r)  dr,  calculation  of  a  transmission  profile  implies  a 
corresponding°extinction  profile.  We  call  this  the  "target"  method  of 
deriving  an  extinction  profile  because  it  requires  the  use  of  a  target 
to  derive  the  oonstant  k/C  value.  Figure  2  presents  an,  example  of  an 
extinction  coefficient  profile  for  a  single  lidar  observation  during  an 
obscurant  test  for  which  the  transmission  history  is  presented  in  Figure 
1.  Figure  3  presents  a  contour  plot  of  obscurant  extinction  (km*1) 
distributions  along  the  observed  path  as  a  function  of  evynt  time. 

If  the  value  of  k/C  remains  oonstant  for  a  given  type  of  obscurant, 
this  value  of  k/C  can  be  used  to  derive  extinction  and  transmission 
values  from  single-ended  lidar  observations  and,  therefore,  vertical  and 
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slant-path  measurements  would  be  possible.  Values  of  k/C  have  been 
calculated  for  each  lidar  shot  for  the  test  with  data  shown  in  Figure  1 
by  using  the  1.06  ya  wavelength  target-derived  transmissions.  In 
addition,  linear  and  best-fit  Gaussian  means  and  standard  deviations  of 
k/C  have  been  calculated.  These  results  are  shown  in  Figures  4  and  5 
where  it  can  be  seen  that  the  two  methods  of  computing  a  mean  k/C  yield 
similar  results,  Flgire  H  shows  that  k/C  can  undergo  both  gradual  and 
rapidly  fluctuating  variability  that  may  be  the  result  of  pulse-to-pulse 
changes  in  lidar  system  performance  (transmitted  energy,  pulca  shape  and 
detector  response)  as  well  as  changes  in  the  optical  properties  of  ,the 
scattering  obscurant.  The  mean  value  of  k/C  was  used  to  recompute  the 
transmission  history  shown  in  Figure  1.  Many  non-physical  transmissions 
(17  percent  of  the  data  record)  resulted  because  equation  1  can  easily 
result  in  negative  transmissions  when  improoer  k/C  values  are  used.'  We 
conclude  that  the  target  method  provides  accurate  lldar-to-target 
transmissions  and  consistent  extinction  coefficient  profiles  when 
multiple  scattering  is  not  Important,  but  that  k/C  must  be  evaluated 
independently  for  each  lidar  observation.  Therefore,  the  method  does  not 
appear  suitable  for  single-ended  vertical  and  slant  path  observations. 


B.  Klett  Method 

The  extraction  of  extinction  coefficients  from  lidar  measurements 
taken  in  an  inhomogeneous  atmosphere,  such  as  through  a  smoke  or  dust 
event,  is  a  difficult  task  in  the  absence  of  a  uniform  reflective 
target.  Early  solutions  to  the  single-scattering  lidar  equation  assumed 
a  backseat ter- to-exti net ion  relationship  and  a  known  value  of  extinction 
at  a  small  distance  from  the  lidar  system  (_*,  5).  However,  as  shown  by 
Klett  (6),  these  forward  integration  solutions  typically  yield  numeri¬ 
cally  unstable  results  unless  unrealistic  solution  parameters  are 
assumed.  This  is  related  to  the  problem  of  evaluating  equation  1  with 
small  uncertainties  in  k/C. 


Klett  overcame  the  stability  problem  by  assuming  a  power  law 
relationship  between  backscatter  and  extinction  and  a  known  value  of 
extinction  near  the  end  of  the  lidar  signature,  leading  to  a  backward 
integration  solution.  In  particular,  the  Klett  solution  for  extinction 
is  as  follows: 


c(R) 


exp  C(S(R)  -  S  )/ k'] 
m 


exp  {[S(r)  -  S  ]/k' }  dr 

Ul 


(2) 


where  S(R)  -  l.itP(R)]t  S  (R  ),  a  -  e(R  ),  and  0  -  const  e  .  The  Klett 

m  oi  m  in 

solution  requires  an  assumed  power  law  exponent  k',  and  knowledge  of  the 

extinction  coefficient  aa  at  a  reference  range  Rm  near  the  end  of  the 

lidar  backscatter  signature.  We  have  applied  the  Klett  solution  with 

several  different  values  of  k*  and  a  to  derive  extinction  profiles  for 

m 

three  different  lidar  shots  made  during  the  obscurant  test  for  which 
transmission  data  are  shown  in  Figure  1. 


For  each  lidar  signature,  extinction  profiles  derived  using  the 
Klett  method  can  be  compared  with  the  extinction  profile  calculated 
using  the  target  method.  Following  the  analysis  of  Sasano  and  Nakane 
(7_),  relative  root  mean  square  errors,  defined  as 


1 

2.2 


Je 


^n  1  Klett (r)  *  "Target^-1  ^ 

r~rl  _ _ 

n  £  aTarget^r* 
r-r^ 


(3) 


where  n  is  the  number  of  range  bins  along  the  path,  were  computed  for  a 
range  of  k*  and  a  values.  These  errors  are  shown  in  Figure  6  as  plots 

ffl 

of  4e  percentage  as  a  function  of  a  and  k'  for  laser  shot  5#6  of  the 

f# 

test  data  shown  in  Figure  1.  The  Klett  extinction  profile  yielding  the 
smallest  ms  difference  from  the  extinction  profile  derived  from  the 
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target  method  is  shown  in  Figure  7.  The  am  and  k'  parameter  used  to 

derive  this  solution  are  indicated  in  Figure  6  by  the  central  x.  The 

near-equal  extinction  profiles  of  the  Klett  and  target  methods  are 

expected,  as  both  solutions  are  based  on  the  single-scattering  lidar 

equation.  Figure  8  presents  comparisons  of  the  target  method  extinction 

profile  with  Klett  method  profiles  for  10  percent  changes  in  k’  and 

factor  of  3  changes  in  a  .  The  values  of  k'  and  a  for  the  four 

m  d 

solutions  shown  in  Figure  8  are  marked  on  the  contour  plot  of  Figure  6 
by  x's  surrounding  the  central  x.  The  changes  in  k'  lead  to  ras  errors 
of  about  30  percent,  and  the  changes  in  o  lead  to  rms  errors  of  nearly 

D 

60  percent.  Figure  9  presents  a  contour  analysis  of  transmission 

evaluated  from  the  Klett  extinction  profile  solutions  as  a  function  of 

the  solution  parameters.  The  transmission  values  for  the  10$  change  in 

k*  and  factor  of  3  change  in  °m  also  are  listed  in  Figures  7  and  8.  A 

10  percent  uncertainty  in  k'  results  in  about  3  percent  error  in 

transmission.  A  factor  of  3  uncertainty  in  a  results  in  about  10 

o 

percent  error  in  transmission.  The  transmission  evaluated  from  the 
target  return  was  19.9  percent,  as  Indicated  by  the  central  x  plotted  in 
Figure  9. 

Figure  10  presents  a  contour  plot  of  rms  differences  between  Klett 
and  target  solutions  for  lidar  shot  number  542.  In  this  case,  the 

transmission  evaluated  by  using  the  target  was  25.6  percent.  Again,  a 

10  percent  change  in  k’  causes  an  rms  error  of  about  60  percent. 

Resulting  extinction  profiles  and  transmission  errors  are  shown  in 
Figure,  11,  with  differences  from  tar get -derived  values  about  the  same  as 
for  laser  shot  546  presented  above. 

Figure  12  presents  a  contour  plot  of  rms  differences  between  Klett 

and  target  solutions  for  lidar  shot  number  550.  In  this  case,,  the 

target  derived  transmission  was  52.3  percent.  The,  10  percent  change  in 

k*  and  faotor  of  3  change  in  o_  give  higher  rms  errors  (nearly  40  and  80 

m 

peroent,  respectively)  than  for  the  lower  transmission  laser  firings 
presented  above.  Resulting  extinction  profiles  and  transmission  errors 
are  shown  in  Figure  13«  Transmission  errors  are  about  4  to  6  percent 
for  10  peroent  ohango  In  k'  and  are  about  17  to  21  peroent  for  a  faotor 
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of  3  change  in  a  .  From  this  result.  It  appears  that  Kiatt  solutions 
m 

for  conditions  of  higher  transmissions  are  more  sensitive  to 

uncertainties  in  k*  and  a  . 

m 


C.  Experimental  Data  Method 

A  number  of  researchers  have  found  a  linear  correlation  between 
backscatter  and  extinction  coefficients  for  typical  particle  size 
distributions  and  visible  or  near-infrared  radiation.  This  implies  a 
linear  relationship  between  optical  depth  and  path- integrated  range- 
corrected  lidar  signal  for  low  density  aerosols.  Because  of 
attenuation,  the  relationship  becomes  nonlinear  for  larger  values  of 
optical  depth  (e.g.,  higher  density  aerosols).  Uthe  (8)  derived  optical 
depth-to-integrated  lidar  signal  relationships  based  on  transmissions 
evaluated  from  target  returns  and  corresponding  lidar  back-scatter 
signatures  for  smoke  and  dust  events.  Using  the  1.06  um  wavelength 
■transmissions  shown  in  Figure  1  and  corresponding  lidar  backscatter 
signatures,  a  plot  of  optical  depth  to  path-integrated  lidar  signal  is 
shown  in  Figure  14.  Each  data  point  is  derived  from  a  single  lidar 
observation  during  the  smoke  event.  A  non-linear  least-squares  curve- 
fitting  method  (9)  was  used  with  the  data  presented  in  Figure  14  to 
evaluate  the  best  fit  a1  and  a2  parameters  of, the  following  expression: 

-a2* 

u(T)  -  a.jT  (1  *  9  )  (4) 

where  T  is  the  path-integrated  lidar  signal  and  u  is  the  optical  depth 

evaluated  from  the  target  return  method.  The  best-fit  relationship  (a1 , 
a2)  and  the  upper  and  lower  uncertainty  limits  as  defined  by  Bevington 

(9)  (a1  ♦  6a1 ,  a2  -  5a,,)  and  (a1  5a^ ,  a2  ♦  6a2)  are  shown  in 

Figure  14. 

For  each  lidar  observation,  an  extinction  profile  can  be  evaluated 
by  first  Integrating  the  lidar  signature  between  the  lidar  location  and 
range  rt  and  then  applying  the  resulting  values  with  the  u(Y) 


relationship  (Equation  4)  to  determine  an  optical  depth  profile  uj. 
Differentiation  yields  the  extinction  profile  o(r).  An  example  of  this 
technique  applied  to  lidar  shots  550,  542,  and  546  (used  for  analysis  of 
the  Klett  method  in  Section  2.2)  is  shown  in  Figure  15.  Solutions  of 
extinction  profiles  derived  by  the  target  method  (Section  2.1)  and  by 
the  experimental  data  method  with  optical  depth-to-integrated  lidar 
signal  expressions  defined  by  (a^i  a2)  ,  (a^  +  5a^ ,  aj  * 

5a£,  and  (a1  -  5a1 ,  a2  +  6a2)  are  shown  in  Figure  15.  These  results 
show  the  experimental  data  method  provides  extinction  profiles  in 
general  agreement  with  the  target  method  for  high  transmissions  but  with 
greater  discrepancy  for  lower  transmissions.  This  may  be  expected  since 
the  target  method  is  based  on  the  single-scattering , lidar  equation, 
while  the  experimental  data  method  includes  effects  of  multiple 
scattering.  Because  of  multiple  scattering,  less  attenuation  correction 
is  applied  at  greater  depths  within  the  smoke  plume  than  with  the  single 
scattering  solution.  For  this  reason,  the  experimental  data  method  may 
be  more  appropriate  for  evaluating  extinction  profiles  through 
relatively  dense  smoke  and  dust  events. 

Ill  CONCLUSIONS 

In  this  paper,  extinction  profiles  have  been  derived  from  lidar 
data  using  three  different  methods.  The  target  method  requires  a  lidar 
target  return  for  each  observation  but  provides  extinction  profiles  that 
when  range  integrated  give  transmission  data  consistent  with  independent 
transmlssometer  observations.  The  Klett  method,  also  based  on  the 
single  scattering  lidar  equation,  can  reproduce  the  extinction  profile 
derived  from  the  target  method  when  appropriate  solution  parameters  are 
used.  However,  expected  uncertainties  in  the  solution  parameters  can 
result  in  unacceptable  errors  in  the  extinction  profile.  The 
experimental  data  method  requires  a  known  relationship  between  optical 
depth  and  path- integrated  lidar  signal.  This  relationship  can  be 
derived  from  lidar  target  and  aerosol  returrts  or  can  be  estimated  from 


9xtinction-to-backscatter  ratios.  This  method  did  not  reproduce  the 
extinction  profiles  of  the  target  or  Klett  methods — probably  because 
this  method  includes  multiple  scattering  effects.  Before  any  one  method 
is  judged  superior,  we  believe  an  experiment  is  needed  that  is  designed 
to  vali iate  lidar  analysis  methods.  Such  an  experiment  will  require  in 
situ  measurement  of  aerosol  optical  and/or  physical  densities  along  the 
propagation  path  of  lidar  observations. 
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FIGURE  1  TRANSMISSION  HISTORY  AND  OPTICAL  DEPTH  RELATIONSHIPS 

FOR  SMOKE  GENERATED  DURING  TRAIL  218  (3RP  ZUNI  ROCKETS) 


FIGURE  2.  RANGE- RESOLVED  VOLUME  EXTINCTION 

COEFFICIENT  (1.06  m)  ALONG  THE  OPTICAL 
PATH  AT  Z  +  63  SECONDS  OF  TRIAL  218 
(TARGET  METHOD) 
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FIGURE  3  CONTOUR  MAP  OF  1.06  pm  WAVELENGTH  EXTINCTION  COEFFICIENT  (km-1) 
DISTRIBUTION  ALONG  THE  LIDAR  OPTICAL  PATH  AS  A  FUNCTION  OF  SMOKE 
EVENT  TIME.  CORRESPONDING  TRANSMISSION  HISTORY  IS  ^LOTTED  ABOVE 
THE  CONTOUR  MAP. 
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FIGURE  7  TARGET  METHOO  AND  BEST-FIT  KLETT  METHOD 
EXTINCTION  COEFFICIENT  PROFILES 
(LASER  SHOT  3461 


EXTINCTION,  o - PER  METER 


FIGURE  9  CONTOUR  PLOT  OF  TRANSMISSIONS  EVALUATED  FROM  TV  E  KLETT  METHOD 
AS  A  FUNCTION  OF  THE  SOLUTION  PARAMETERS  and  K'.  THE  X'l 
CORRESPOND  TO  THOSE  OF  FIGURE  S. 
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FIGURE  13  EXTINCTION  COEFFICIENT  PROFILES  FOR 
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Appendix  B 


"ERROR  ANALYSIS  OF  LIDAR  SOLUTION  TECHNIQUES  FOR 
RANGE-RESOLVED  EXTINCTION  COEFFICIENTS  BASED  ON  OBSERVATIONAL  DATA" 

Presented  at  SMOKE/OBSCURANTS  SYMPOSIUM  IX. 
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ERROR  ANALYSIS  OF  LIDAR  SOLUTION  TECHNIQUES  FOR  RANGE-RESOLVED  EXTINCTION 
COEFFICIENTS  BASED  ON  OBSERVATIONAL  DATA* 

E.  E.  Uthe  and  J.  M.  Livingston 
SRI  Intarnatlonal 
Menlo  Park,  CA  9*025,  USA 


ABSTRACT 


The  SRI  four-wavelength  (0.53,  1.06,  3.8  and  10.6  imO  Ildar  system  wss  used  during  the  Siaoke 
Week  XI /SNOW-TWO  field  experiments  to  observe  range-resolved  backscatter  and  target-reflected  laser 
energy  along  a  horizontal  path  intersected  by  test  obscurant  aerosols. 

Target  reflected  laser  energy  provided  an  excellent  measure  of  time-dependent  path-integrated 
extinction  (e.g.,  transmission).  These  values  allow  the  range-resolved  backscatter  to  be  analyzed 
In  terms  of  range-resolved  extinction  when  .-ingle-scattering  is  assumed.  This  two-ended  (lidsr  and 
target)  method  provides  extinction  ost-tf lclents  along  the  optical  path  that  are  consistent  with 
Independent  transmlssometer  otrar vat ions.  The  values  are  the  most  accurate  that  can  be  derived  from 
the  Ildar  technique  when  multiple  scattering  Is  limited. 

Extinction  profiles  evaluated  from  the  lldar/target  method  can  be  utilized  to  Investigate  the 
limitations  and  accuracies  of  algorithms  which  are  used  to  derive  extinction  data  rrom  Ildar 
backscatter  signatures  when  no  target  return  Is  available.  Single-ended  Ildar  methods  are  needed 
for  vertical  and  slant-path  measurements.  Two  single-ended  techniques  were  investigated.  The  first 
is  the  well-known  Klett  modification  of  the  analytic  solution  to  the  single-scattering  Ildar 
equation.  The  Klett  technique  can  reproduce  the  extinction  profile  derived  from  the  lldar/target 
method  as  expected  since  both  are  derived  from  the  single-scattering  lldar-equatlon.  However,  this 
requires  relatively  accurate  Information  on  a  boundary  value  of  extinction  and  the  backsoatter-to- 
extlnctlon  relationship.  Contour  plots  showing  deviation  of  the  Klett  extinction  profile  from  the 
lldar/target  extinction  profile  as  a  function  of  the  solution  parameters  are  presented. 

A  second  method  of  deriving  extinction  profiles  from  single-ended  Ildar  measurements  Is  based 
on  an  observed  relationship  of  optical  depth  to  path-integrated  Ildar  slpial.  £:ich  a  relationship 
derived  from  lldar/target  data  lnoludes  multiple  scattering  effects  on  the  extlncci^n  evaluation. 
The  method  does  not  reproduce  the  extinction  profile  derived  from  the  single-scattering  lldar/target 
and  Klett  methods  —  Indicating  the  importance  of  multiple  scattering  on  the  Ildar  backscatter 
signature. 

1.  BACKGROUND 


Description  of  the  SRI  van-mounted  four-wavelength  (0.53,  1.06,  3.8  and  10.6  w»)  Ildar  system 
and  its  operation  during  the  Smoke  Week  VI /SNOW-TWO  field  program  was  presented  at  the 
Smoke/Obsourants  Symposium  VIII  (DeLateur  et.al.,  198*).  Basically,  laser  pulses  at  the  four 
wavelengths  were  transmitted  along  a  horizontal  path . terminated  at  a  distance  of  600  a  from  the 
Ildar  by  a  target  of  uniform  reflective  properties  over  the  illuminated  area.  The  laser  pulses  at 
each  wavelength  were  emitted  sequentially  within  a  60  pa  interval  to  reduce  effects  of  time-  and 
space-varying  aerosol  concentrations  on  each  four-wavelength  pulse  set.  About  60  pulse  sets  were 
transmitted  for  each  minute  of  operation.  Collected  energy  backvcattered  from  atmospheric  aeroso’w 
und  energy  reflected  from  the  solid  target  was  collected  within  a  telescope  and  directed  to 
appropriate  detectors  using  dlchrolo  filters.  After  logarithmic  amplification,  backscatter 
signatures  were  digitised  by  an  8-bit  100-MHa  transient  recorder  and  resulting  digital  records  were 
sttred  on  nine-track  magnetic  tape  for  subsequent  analysis. 


•This  work  was  supported  by  the  U.S.  Ai-my .  Research  Office,  Geosciences 
Division,  under  Contract  DAAG29-82-K-0191. 
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2. 


DATA  ANALYSIS  METHODS  AND  RESULTS 


2.1  Target  Method 

Target  returns  observed  by  a  Ildar  system  can  be  interpreted  in  terms  of  transmission  across 
the  lldar-to-target  path.  Previous  results  have  shown  that  l.dar  derived  transmissions  agree  with 
independent  data  records  collected  with  conventional  two-ended  transmissoeieters  (Uthe  et  al., 
1983).  Therefore,  the  four-wavelength  Ildar  target  returns  oollected  during  an  obscurant  event 
pro.  .de  the  time  history  of  transmission  at  aach  wavelength  and  provide  data  to  evaluate  optical 
depth  u«-ln  T  relationships  between  wavelengths,  as  shown  by  the  example  presented  in  figure  1. 


0  1  23  01  2  3  01234 

OPTICAL  DEPTH  <»T  INOICATEO  WAVELENGTH  (X)  —  Uv 

FIGURE  1  TRANSMISSION  HISTORY  ANO  OPTICAL  OEPTH  RELATIONSHIPS 

FOR  SMOKE  GENERATED  DURING  TRAIL  218  (3RP  ZUNI  ROCKETS) 


Fernald  et  al.  (1972),  among  others,  have  shown  that  If  the  ratio  of  the  extinction  coefficient 
to  the  backscatter  coefficient,  k-«(r)/S(r),  remains  constant  over  the  range  Interval  A r  being 
Investigated,  then  the  two-way  transmittance,  T^r),  follows  directly  from  differentiates  of  the 
basic  Ildar  equation  and  Is  given  by 


T^(R) 


-  ^  J  P(r)dr 
C  0 


(1) 
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In  Uilf  eouatisn,  P(r)  ia  the  signal  whose  amplitude  is  proportional  to  the  power  received  from  a 
scattering  voluoe  at  ranga  r  and  is  oo rrectad  for  the  geometric  range-squared  affect  and  transacted 
*”«r|y  variations.  C  is  a  Ildar  calibration  oonstant. 

'ha  constant  k/C  can  be  determined  by  combining  the  known  two-way  transaittance  between  the 
laser  and  the  target  with  the  Ildar  return  signature  integrated  over  the  saae  range.  The  value  of 
tlv»  const/.'  ■  '  Z  .a  ,  en  substituted  back  into  (1)  and  ooablned  with  Integrals  of  the  received 
backscatter  alc-’Als  to  intermediate  ranges  to  yield  a  horizontal  transaisslon  profile.  Since 
transalssi"  ^s  directly  related  to  optical  depth  through  T2!*)  -  e-2u^R\  where  optical 
depth  u(Rj  >  f  s(r)  dr,  calculation  of  a  transaisslon  profile  laplles  a  corresponding  extinction 
profile,  We  3ll  this  the  "target*  sethod  of  deriving  an  extinction  profile  because  it  requires  the 
use  of  i  target  to  derive  a  constant  k/C  value.  figure  2  presents  an  example  of  an  extinction 
coefficient  profile  for  a  single  lldrr  observation  during  an  obscurant  test  for  which  a  transaisslon 
history  is  presented  in  figure  1. 

If  the  value  of  k/C  remains  constant  for 
a. given  type  of  obscurant,  this  value  of  k/C 
can  be  used  to  derive  extinction  and 
transaisslon  values  froa  single-ended  lidar 
observations  and,  therefore,  vertical  and 
slant-path  aeastreaents  would  be  possible. 

Values  of  k/C  have  been  calculated  for  each 
Ildar  shot  for  the  test  with  data  shown  in 
figure  t  by  using  the  1  06  as  wavelength 
target-derived  transmissions.  In  addition, 
linear  and  beet-fit  Gaussian  means  and 
standard  deviations  of  k/C  have  been 
calculated.  These  results  are  shown  in 
figures  3  and  *  where  it  can  be  seen  that  the 
two  aethods  of  oomputlng  a  mean  k/C  yield 
similar  results,  figure  3  shows  that  k/C  can 
undergo  both  gradual  and  rapidly  fluctuating  variability  that  may  be  the  result  of  pulse-to- pulse 
changes  in  Ildar  systea  performance  (tranaaitted  energy,  pulse  shape  and  detector  response)  i  well 
aa  manges  in  the  optical  properties  of  the  scattering  obscurant.  The  bean  value  of  k/C  was  used  to 
recompute  the  transmission  history  ahowi  in  figire  t.  Many  non-physical  transmissions  (IT  percent 
of  the  data  reoord)  resulted  because  equation  t  can  easily  result  in  negative  transmissions  when 
Improper  k/C  values  are  u  *1.  We  oonclude  that  the  target  sethod  provides  accurate  lldar-to-target 
transmissions  and  consistent  extinction  coefficient  profiles  when  multiple  test taring  is  not  impor¬ 
tant,  but  that  k/C  must  be  evaluated  independently  for  seen  lidar  observation.  Therefore,  the  sethod 
does  not  appear  suitable  for  single-ended  vervloal  and  slant  patn  obaervitiona. 

2.2  Klstt  Method 

The  extraetlyn  or  extinction  coefficients  from  lidar  as aaur amenta  taken  in  an  inhomogeneous 
atmosphere,  such  as  through  a  smoke,  or  dust  "vent,  U  a  difficult  task  tn  the  abeenoe  of  a  uniform 
reflective  target.  Carly  solutions  to  the  single-eeattering  lidar  equation  assumed  a  backseat t #r- 

i 
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FIGURE  2.  RANGE-  RESOLVED  VOLUME  EXTINCTION 

COEFFICIENT  (l  OSiewl  ALONG  THE  OPTICAL 
PATH  AT  Z  ♦«  SCCONOS  OF  TRIAL  >.* 
(TARGET  METHOOI 


A-2S 


to-extlnction  relationship  and  a  known  value  of  extinction  at  a  small  distance  from  the  Ildar  system 
(V iezee  et  al.,  1969 ;  Johnson  and  Uthe,  19T1).  However,  an  shown  by  lUett  (19611.  these  forward 
Integration  solutions  typically  yield  numerically  unstable  results  unless  unrealistic  solution 
parameters  are  assumed.  This  la  slaillar  to  the  pro  ties  of  evaluating  equation  i  with  small 
uncertainties  in  k/C. 


SHOT  MUMMER 


FtOUhf  J.  K/C  MISTOhV  SOS  MOKt  QENf RATIO  DURING  TRIAL  Hi 


0  02  04  0-S  OS  IS 

UOAR  CONST  AMT.  k/CalO"* 

fhguM  4  otrrnieuTxx  or  kk  vsluU  son 
MSOKI  TRIAL  JH . 

Klett  overcame  the  stability  problem  by  assuming  a  power  law  relationship  between  backscatter 
and  extinction  and  a  Mown  value  of  extinction  near  the  end  of  the  Ildar  signature,  leading  to  a 
backward  integration  solution.  In  particular,  the  Klett  solution  for  extinction  is  aa  foil owe ■ 
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*xp  C(S(R>  -  SVk’3 

•  (R)  -  - -x- - - - ,  (2) 

•*’  ♦  £,  J*  .XP  {CS(D.  -  S  Vk-|  dr 

*  « 

k' 

where  S(R)  •  ln[P(R)],  *  SCR,,),  »m  -  atR^) ,  and  t  •  coast  e  .  The  Klett  solution  requires  an 

mi— rl  powar  law  axponant  k* ,  and  knowladg*  of  the  extinction  coefficient  at  the  reference  range 
He  have  applied  the  Klett  solution  with  several  different  values  of  k'  and  a  to  derive 
extinction  profile*  for  three  different  Ildar  shots  aade  during  the  obscurant  test  for  which 
transalsalon  data  Is  shown  In  Figure  1 . 

For  each  Ildar  signature,  extinction  profiles  derived  using  the  Klett  set nod  can  be  compared 
with  the  extinction  profile  calculated  using  the  target  eethod.  Following  the  analysis  of  Sasano 
and  Kakane  (1984),  relative  root  aean  square  errors,  defined  as 


1 

2l2 


M  • 


H?  *  Klett(r)  ■  ‘target (p)^ 
1 


1  trV 


(r) 


(3) 


where  n  Is  the  nuaber  of  rang*  bins  along  th*  path,  were  ooaputed  for  a  rang*  of  k'  and  oa  values. 


F0WIR  UMS  IXFONgNT,  K’ 

FIOUM  |  CONTOUR  nor  OF  kutt  MCTMOO  MWg  (SSOttt  At  A  function  of  tm« 
•Mjution  fanaantim  r  ano  *.  (lams  mot  MSI 
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These  errors  art  shown  in  Figure  5  as  plots  of  is  percentage  as  a  function  of  ea  and  k’  for  laser 
snot  5*6  of  the  teat  data  shown  in  Figure  1.  The  Klett  extinction  profile  yielding  the  smallest  ms 
difference  froa  the  extinction  profile  derived  free  the  target  aethod  is  shown  in  Figure  6.  The 
near-equal  extinction  profiles  of  the  Klett  and  target  nethods  are  expected  as  both  solutions  are 


figure  •  target  methoo  ano  kst-fit  klett  methoo 

EXTINCTION  COEFFICIENT  PROFILES 
(LASER  SHOT  MS) 

uncertainty  in  transmission.  A  factor  of  3 
uncertainty  in  sa  results  in  about  10  percent 
uncertainty  of  transmission.  The  trans¬ 
mission  evaluated  from  the  target  return  was 
19.9  percent. 

Figure  8  presents  a  oontour  plot  of  ras 
differences  between  Klett  and  , target  solu¬ 
tions  for  Ildar  ahot  number  5*2.  In  this 
oaae.  the  transmission  evaluated  by  using  the 
target  was  25.6  percent.  Again,  a  10  percent 
ohange  in  k'  causes  an  ran  error  of  about  60 
percent.  Resulting  extinction  profiles  and 
transmission  uncertainties  are  shown  in 
Figure  9  and  dlfferenoee  from  target-derived 
values  are  about  the  sane  aa  for  laser  ahot 
5*6  presented  above. 

Figure  IQ  presents  a  oontour  plot  of  me 
dlfferenoee  between  Klett  and  target 
solutions  for  Ildar  shot  number.  $50.  In  this 
oeee,  the  target  derived  transmission  was 
52.3  percent.  The  10  percent  ohange  in  k' 
and  factor  of  3  ohange  in  «a  give  higher  me 
errors  (nearly  *0  and  80  percent 
respectively)  than  for  the  lower  transmission 


based  on  the  single-scattering  Ildar  aqua¬ 
tion.  Figure  T  presents  comparisons  of  the 
target  aethod  extinction  profile  with  Klatt 
method  profiles  for  10  percent  changes  in  k* 
and  factor  of  3  changes  In  a  The  values  of 
k*  and  «a  for  the  four  solutions  shown  in 
Figure  ?  are  aarked  on  the  contour  plot  of 
Figure  5  by  x's.  The  changes  in  k'  lead  to 
ras  errors  of  about  30  percent  and  the 
changes  In  «a  lead  to  ras  errors  of  nearly  60 
peresnt.  The  resulting  transmission  errors 
are  shown  in  Figure  T.  A  10  percent  uncer¬ 
tainty  in  k*  results  in  about  3  percent 


(LASiR  SHOT  MSI. 
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POWER  LAW  EXPONENT,  *• 

FIGURE  ■  CONTOUR  PLOT  OP  KLETT  MCTMOO  RMS  ERRORS  AS  A  FUNCTION  OP  THE 
SOLUTION  PARAMETERS  K’  AMO  «„  (LASER  WOT  M2) 


Lamp  firings  presented  above.  Resulting  extlnetlon  profiles  and  tranaslsalon  uncertainties  ere 
•ham  la  rigure  11.  Trenselselon  differences  irt  about  *  to  6  percent  for  10  paroant  ohange  in  k’ 
and  are  about  IT  to  21  paroant  for  a  factor  of  3  oftanca  In  if.  Pro*  this  result,  it  appears  that 
Klatt  aolutiona  for  oondltlona  of  higher  traneslssions  ara  aora  aanaltlva  to  uncertainties  in  k* 
and  e_. 


2.3  . Ixperlsental  Data  Method 

A  nunbar  of  raaaarchara  haaa  found  a  linear  correlation  between  back acat tar  and  axtlnotlon 
coefficient*  for  typloai  particle  alM  dlatributtona  and  visible  or  neor-infrved  radiation.  Thla 
liplUi  a  linear  ralationahlp  between  optloal  depth  and  path* integrated  range-oorreoted  Ildar  signal 
^or  density  aerosols.  tithe  0960  derived  optloal  depth -to- integrated  Hoar  signal 
relationships  based  on  tranaalsslona  evaluated  fro*  target  returns  and  corresponding  Ildar 
baoksoatter  signatures  for  awoke  and  duet  events.  Using  the  1 .06  wavelength  tranaalsslona  shown 
In,  Figure  1  and  oorrespondlng  Ildar  baoksoatter  signatures,  a  plot  or  optloal  depth  to  path 
integrated  Ildar  signal  la  shorn  In  Figure  12.  Each  data  point  la  derived  fro*  a  single  Ildar 
observation,  a  non- linear  least-squares  ourve-flttlng  aethod  (levington.  1969)  was  used  with  the 
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-  data  presented  Xn  Figure  12  to  evaluate  the  bast  fit  a,  and  a2  paraaetars  of  the  following 
expression: 


u(T) 


*2* 

a.T  (1  ♦  •  e  ) 


(*> 


where  T  Is  ths  path- Integrated  Ildar  signal  and  u  Is  tha  optical  dapth  evaluated  from  the  target 
return  aethod.  The  best-fit  relationship  (a,,  a2)  and  the  upper  and  lower  uncertainty  Halts  as 
defined  by  Bsvlngton  (a,  ♦  4a, .  a2  -  «a2)  and  (a,  -  4a, ,  a2  ♦  4a2)  are  shown  In  Figure  12. 

For  each  Ildar  observation,  an  extinction  profile  can  be  evaluated  by  first  Integrating  the 
Ildar  signature  between  the  Ildar  location  and  range  r,  and  then  applying  the  resulting  v,  values 
with  the  u(T)  relationship  (Equation  *)  to  determine  an  optical  depth  profile  Uj.  Differentiation 

yields  the  extinction  profile  o(r).  An 
example  of  this  technique  applied  to  Ildar 
shots  550,  542,  and  5*6  (used  for  analysis  of 
the  Klett  aethod  In  Section  2.2)  Is  shown  in 
Figure  13.  Solutions  of  extinction  profiles 
derived  by  the  target  aethod  (Section  2.1) 
and  by  the  experlaental  data  nethad  with  op¬ 
tical  depth-to-integrated  Ildar  signal 
expressions  defined  by  (a,,  a2),  (a,  ♦  4a, , 
a2  -  4a2).,  and  (a,  -  4a, ,  aj,  *  4a2)  are  shown 
In  Figure  13.  Theea  results  show  the  experi- 
aental  data  aethod  provides  extinction  pro- 
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files  In  general  agreement  with  the  target 
aethod  for  high  transmissions  but  with 
greater  dlsorspency  for  lower  tranaals- 
alona.  This  nay  be  expected  since  the  target 
aethod  Is  based  or  the  single-scattering 
Ildar  equation  while  the  experlaental  data 
aethod  lnoludea  effects  or  multiple  scat¬ 
tering.  For  this  reason,  the  experlaental 
data  aethod  aay  be  aore  appropriate  for 
evaluating  extinction  profiles  through 
relatively  dense  aaoke  or  dust  events. 


Ftouag  g  extinction  coefficient  frofh.ES  for 
INCREMENTED  VALUES  OF  K*  ANO  tm 
(LASER  SHOT  M2 1 


3.  CONCLUSIONS 


In  this  paper,  extlnotlon  profiles  have  been  derived  froa  Ildar  data  using  three  different 
aethods.  The  target  aethod  requires  a  Ildar  target  return  for  each  observation  but  provides 
extlnotlon  profiles  that  when  range  integrated  gives  transmission  data  consistent  with  Independent 
tranaaiaaoaeter  observations.  Ths  Klett  aethod,  also  baaed  on  the  tingle  test taring  ltdar  equation, 
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can  reproduce  the  extinction  profile  derived  free  the  target  eethod  when  appropriate  solution 
parameters  are  used.  However,  expected  uncertainties  In  the  solution  parameters  can  result  in 
unacceptable  errors  In  the  extinction  profile.  .  The  experimental  data  method  requires  a  known 
relationship  between  optical  depth  derived  from  Ildar  target  returns  and  path- integrated  Ildar 
signal  derived  from  aerosol  scattering.  Extinctlon-to-backscatter  ratios  can  provide  estimates  of 
the  needed  relationships.  This  method  did  not  reproduce  the  extinction  profiles  of  the  target  or 
Klett  methods — probably  because  this  method  Includes  multiple  scattering  effects.  Before  any  one 
method  is  Judged  superior,  we  believe  an  experiment  Is  needed  that  Is  designed  to  validate  Ildar 
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FOOR-WAVELENGTH  LIDA*  MEASUREMENTS 
non  SMOKE  WEEK  71/SHOW-TWO 

S.A.  DeLateur,  H.B.  Nielsen, 

E.E.  Ctha  and  J.M.  Livingston 
SEX  International 
Manlo  Park,  CA  9*025 

ABSTRACT  (U) 

(0)  SIX  Intaraa clonal 'a  four -wavelength  Ildar  system  was  ua ad  to  aaka  backacattar  and  trans¬ 
mission  observations  at  Xaaar  wavelength*  of  0.53,  1.96,  3.8,  and  10.6  v*  daring  eh*  Smoke  Waak  VI/ 

SHOW-TWO  {laid  r-ograa.  Laaar  pula a a  wore  traaaaietad  along  a  horizontal  path  tarainacad  at  a  dlatanca 
of  600  a  from  tba  Ildar  by  a  targat  with  uniform  reflective  propartlaa  over  Cba  lllualnatad  araa. 

Targac  returns  ara  lntarpratad  In  terms  of  claa-depeadant  laaar  energy  cranaalaalon  ovar  tba  path,  and 
provide  tha  aaana  to  lntarprat  lidar  backacattar  ilgnaturaa  In  taraa  of  ranga-dapandant  axtlnctlon 
eoafflclanta.  Although  tha  data  aaalyala  program  la  in  ita  early  ataga,  an  example  of  tha-  technique  la 
praaanead.  Soma  raaulta  obtalnad  from  tha  SWOW-OBX-I  taata  ara  alao  lncludad. 


1.  (U)  ISTRODUCTIOB 


(U)  Tha  Saoka  Waak  VI /SHOW- TWO  program  «u  conducted  by  tha  U.S.  Army  Saoka/Obacuranta  PM  and 
CRKEL  during  tha  1983-8*  winter  at  Grayling,  Michigan.  Varloua  optical  aanaora  warn  used  to  sake  propa¬ 
gation  maaauraaanta  throrgh  both  natural  and  nan -mad  a  obacuranta.  With  funding  from  tha  U.S.  Army 
kaaaarch  Of  flea,  SKI  Iaearnaclonal  pavtldpatad  by  naklng  lidar  aaaauramanta  alaultanaoualy  at  wave- 
langtha  of  0.53,  1.06,  3.8,  and  10.6  ua. 


(U)  Tha  objectlvaa  of  thla  program  vara  to 

a  Evaluate  tranaalaalona  for  natural  and  generated  obacuranta  at  DoD  laaar  wavelengtha 
■  a  Evaluate  nultlpla-acaetarlng  affects  with  nultl-vavalaagth  lidar  obaervaclona 
a  lntarprat  nultl-wavalangth  backacattar  slgnsturea  In  terms  of  aaroaol  optical,  and 
physical  danaltlea,  and  particle  characteristics,  along  tha  propagation  paths 
a  Derive  statistical  aodals  describing  tha  special  and  temporal  variability  of  obscurants 
along  propagation  paths. 


2.  (U)  MULTIWAVELEMGTH  LUJAK  STST2M 


(U)  Tha  mulclwavelaagth  lidar  syatsm  Is  Installed  within  a  6-mater  long  van  to  facilitate  opera¬ 
tions  at  remote  sitae  (Figure  1).  Laser  pulses  ara  transmitted  either  vertically  or  horizontally  at 
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FIGURE  1  FOUR-WAVELENGTH  LIOAR  SYSTEM 
IN  8-METER- LONG  VAN 


wavelengths  of  0.53,  1.06,  3.8  and  10.6  wa  by  using 
Md:YAG,  DF,  and  CO^  laaara.  Tha  cranaalttar  optics 
direct  tha  laser  pulses  to  a  co— on  path  coaxial 
with  a  12-in  Hawtonlan  tala sc ope,  a a  shown  in 
Figure  2.  A  series  of  laser  pulses  at  each  of  the 
four  wavelengths  are  transmitted  sequentially  with¬ 
in  a  60  us  Interval.  Therefore,  effects  of  tine- 
end  space-varying  aerosol  distributions  on  tha 
aultl-wavelength  Ildar  measurements  are  nlnlalzed. 


(U)  Collected  energy  baefcseattered  f row  atmo¬ 
spheric  targets,  and  energy  reflected  from  solid 
targets.  Is  collected  within  the  telescope  and  di¬ 
rected  to  appropriate  detectors  using  dlchrolc 
filters.  Detector  outputs  are  logarlthalcally  am¬ 
plified  and  digitised  by  8-bit  100- MHo  and  12-blt 
10-MHs  transient  recorders.  Transient  recorder  out¬ 
puts  are  displayed  on  an  oscilloscope,  and  are 
written  on  magnetic  tape  under  control  of  an  LSI-11 
microcomputer  system. 

(TJ)  The  four-wavelength  system  was  first  used 
during  tha  S NOV-O NE-B  tests  conducted  at  Grayling, 
Michigan,  during  November  and  December  1983.  Data 
collected  showed  that  greater  receiver  dynamic  range 
and  shorter  recovery  time  from  large  amplitude 


FIGURE  2  BLOCK  DIAGRAM  OF  THE  SRI  FOUR-WAVELENGTH  LIOAR 
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signals  war*  needed  to  evaluate  dsns*  obscurant  distributions  along  tha  propagation  path.  A  significant 
part  of  our  prograa  was  diractad  to  iaprovaasnt  of  tha,  lidar  racaivara  for  application  to  tha  Saok*  Hack 
VI/SNOW- TWO  experiments.  Saw  extended-range  logarithaic  amplifiers  war*  i  *railad  on  the  0.S3  urn  and 
1.06  ua  wavelength  channels.  Tha  now  amplifiers  aztandad  tha  dynaalc  rang*  froa  four  to  six  order*  of 
aagnituda  (i.a.,  froa  40  to  60  dB).  Bacauaa  tha  dat«ctor  foi  tha  3.8  and  10.6  ua  wavalasgth  channel* 
was  tha  Halting  component,  a  tlaa-gata  circuit  was  addad  to  tha  daeactor  pra-aaplifier  so  that  danaa 
aarosol  returns  near  tha  lidar  could  be  electrically  attenuated ,  while  returns  at  greater  distances 
could  be  fully  amplifier*. 

CO)  In  addition  to  the  receiver  aodiflcatlona,  tha  coaputer-con  troll  ad  data  digitisation  and  coe- 
ponants  in  the  computer  systaa  ware  updated.  New  8-bit  100-MHx  and  L2-bit  10-MHs  transient  recorders 
vara  installed  within  a  CAMAC  chaaala  along  with  tha  LSI-11  coaputar.  Softvara  was  stored  on  5-1/4  inch 
floppy  disk  units  which  vara  also  ins tall ad  within  tha  CAMAC  chassis.  Therefore,  tha  eye ter  was  sub¬ 
stantially  reduced  in  sisa  while  increasing  capabilities  and  lap  rowing  rejection  of  BFI  noise  fron  laser 
firing*.  Tha  Saok*  Weak  VT/SSOW-TMO  prograa  waa  tha  first  teat  of  thaaa  new  lidar  aodiflcatlona. 

3.  (0)  DATA  COLLECTION  PBOGBAH 

(O)  Tha  faur-vavelsngth  lidar  wan  waa  transported  to  Grayling.  Michigan,  and  positioned  alongside 
transaisaoaeter  eye  teas  oparatad  by  tha  U.S.  Army  A  tao  spheric  Science  Laboratory  (ASL)  and  tha  U.S. 

Mary  (I0tL).  A  solid  lidar  earget  was  installed  about  600  n  distant  froa  tha  lidar,  and  near  tha  psth- 
tendnatlou  hardware  for  tha  tranaaiasaaatar  sensors.  Therefore,  tha  Ildar-target  path  was  nearly 
parallel  and  saparatad  by  only  a  few  asters  froa  tha  tranaalssoastsr  paths. 

(O)  Lidar  observational  periods  wars  restricted  in  duration  bacauaa  of  laser  cya-safaty  considera¬ 
tions.  Operations  ware  conducted  only  at  flaws  whan  approval  was  grant ad  by  tha  field  coordinator. 
Bacauaa  of  these  restriction*,  only  Halted  lidni  calibration  data  ware  obtained.  Tha  lidar  normally 
oparatad  about  five  minutes  bafora  tha  start  of  each  obscurant  tast  until  flve-to-ren  ainutss  sf tar  tha 
tast  start.  Unlike  lidar  operations  on  tha  SSOW-OW  I  tssts,  laser  transmissions  at  all  four  wavelengths 
ware  allowed.  Tha  lidar  operational  class  sad  obscurants  observed  are  given  in  Table  I. 
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tabu  i.  (j)  lxdar  operaitohs 

Sack*  U««k  VI/SNOW-TUO 


TRIAL 

DATE 

TIME 

TAPE 

MATERIAL  TTPB 

101 

1/10/84 

10:42-10:50 

6 

Diesel  Fuel  Generator 

102 

1/11/84 

11:07-11:16 

10 

fog  Oil  Generator 

202 

1/U/S4 

.  11:*2-U:55 

10 

fog  Oil  Generator 

114 

1/11/L4 

13:58-14:06 

11 

RP  froa  L8A3 

214 

1/11/84 

14:38-14:46 

11 

RP  froa  L6A3 

203 

1/12/84 

8:58-9:12 

12 

IR2  Generator 

ua 

1/12/84 

U:33-U:45 

12 

RP  Zunl  Rocker 

'  210 

1/12/84 

12:00-12:2* 

13 

PUP  2 uni  Rocket 

218 

1/12/84 

14: 38-13  :U 

13 

(3)  RP  2unl  Rockets 

107 

1/12/84 

15:43-13:36 

14 

foreign  8C  Pots 

105 

1/12/84 

16:56-17:02 

13 

Experimental  "C"  Generator 

122 

1/13/44 

8:38-19:16 

16 

fog  Oil  Generator  in  Snow 

222 

1/13/84 

9:28-9:36 

16 

fog  Oil  Generator  in  Snow 

10* 

1/13/84 

13:03-13: U 

17 

1X2  ♦  Fog  Oil  Generator 

204 

1/13/84 

13: 21-13: *0 

17 

IR2  *  fog  Oil  Generator 

U6 

1/13/84 

14:24-14:31 

18 

Slaulatad  135aa  HE 

310 

1/13/84 

13:03-13:22 

18 

(3)  PUP  Zusl  Rockets 

418 

1/13/8* 

13:33-16:08 

19 

(3)  RP  Zunl  Rockets 

111 

1/13/8* 

16:49-17:02 

19 

UP  LOST 

216 

1/14/84 

9:06-9:10 

20 

SlauUtad  133m  HE 

213 

1/14/84 

9:37-10:01 

20 

1X2  XM76 

123 

l/li/84 

8:39-9:1* 

a 

Experimental  "C“  Generator 

128 

1/16/8* 

12:30-12:47 

22 

foreign  Pots 

302 

1/16/8* 

13:45-13:53 

23 

fog  Oil  Generator 

318 

1/16/84 

16:30-16:42 

23 

SlauUted  135m  HE 

4.  (B)  DATA  RICOSDJ 

(0)'  A  four  mnlitk  koduuttir  dlfltii  record  U  tbon  U  figure  3.  The  (ruilMC  tteordtr  la 
Initiated  just  prior  to  the  flrac  User  firing,  and  digitises  continuously  sc  10  ns  (1.3a)  lacarrals 
for  above  AO  us.  Tha  clas  period  beeveaa  tbs  0.S3  m*  ’.0#  uo  wavelength  4s ts  Is  determined  by  ebo 
•spsMCtM  tlaa  between  too  Usor  pulses  dartre*  froa  the  saao  flask  lsap  aaclcaeioe  of  ebo  Md:YAG  Usor 
A  relatively  loot  period  was  used  bo  trees  the  3.0  ua  sad  10.6  ua  User  firings  as  ebo  seas  detector  is 
need,  aad  aaple  tlaa  is  needed  for  detector  stabilisation  after  observing  the  Urge  aaplleude  3.1  ua 
target  retun.  figure  3  shoes  thee  the  1.3  a  (100-Mb)  raeolucloa  of  ebo  8-bit  digitiser  provides  valid 
target  aaplltude  tlfnaU  sc  all  wavelengths.  The  15-e  resolution  (KHOs)  of  the  12-kit  digitiser  did 
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•*  aal  table  cart  at  aapllcudo  4iu  at  0.53  ua  aad  1.06  m.  Accordingly,  cha  8-blt  10O-MU  data 

tacndiiii  ara  balag  aaad  oa  cha  data  aaalyaia  pragma. 

(0)  Tha  foar-aaualaagth  bockacattor  algaatara  (tor ad  alt him  Cha  craaalaat  racordar  taor r  (FI jura 
3)  caa  ba  aaad  to  laaarata  individual  «rral«a*th  bockacattor  algaaturaa  oa  cha  dlaplay  oacllloacopa  (»aa 
Flfara  1).  figure  A  praaabta  aocUlooeopo  diaplara  for  claar  alt  coadltloaa  oa  13  January.  Tha  Initial 
sl^al  rlaa  la  a  raaalt  o t  tha  caawor gaaco  at  tha  traaaalttad  galea  with  tha  roaalvar  field  of  rto*. 
tba  claar- air  raago  raaaload  al«aal  falla  off  aa  laueroa  raaga-ogaarad.  Tha  target  ratura  la  600. a  tram 
tba  Udar  pa  alt  too.  Tba  ratata  at  abaac  300  a  raaalta  tram  a  gar  alta  aaad  to  a  of  part  tha  aataorvloglcal 
tDW,t*  W  •**»  »aa  laatallad  af tar  Ildar-target  poeltloalag,  aad  June  bafora  tba  Saab a  Waak  71 
aartaa  of  too  to.  Wire  aad  targat  rataraa  grande  laforaatloa  oa  palaa  abapa  at  oacb  laaar  uavelaagtb. 

Tha  abort eat  laaar  palaa  daracloa  la  at  0.33  ua  follovad  bp  tba  1.00,  10.6  aad  3.6  ua  uavelaagtb  palaaa. 

(U)  Sacaaaa  tha  elaer-air  nattering  daaraaaaa  with  laaraaalag  uaaalaagtb  ufclle  tba  targat  ratuna 
raaaU  relatively  eoaacaat  with  uavelaagtb,  tba  target-te-eiaer  air  ratio  laaroaaaa  with  wavelength  aa 
ahaaa  aa  tba  0.33  aad  1.06  «  uaaalaagtb  elgaaturoe.  Claar  air  aaattarlag  at  3.6  aad  10.6  ua  la  obaorrad 
oaly  at, abort  raagaa.  Tbarafaro.  obaarratlaa  of  claar  air  aaattarlag  aad  targat  rataraa  reguirae  a 
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CLEAR  AIR  13  JANUARY  1384 

8  lit  -  100  Wz 


100  n/oiv  .532m**  100  *V»iv  1.06m** 


100  n/div  3.8  m**  100  n/oiv  10.6  m** 

FKJUkf  4  LIOAN  •ACXSCATTtR  HQMATUMg  AT  OJ1, 3.6  and  104  an*  FOA  CLgAN  Ah. 
CONOiTKX* 


grantor  nctlw  dynamic  w<i  at  laager  aaial  i'r>«  «•  MfoimM  la  geetloe  2,  tlM  dytuanic  rang*  of 

the  0.33  aa A  l.M  m  oareleagtk  racalwr*  mm  lowtuM  hy  add lag  wraNid  Ayaaalc  raage  (60di)  loga- 
rlUate  aapll flare.  A  fmallflM  gat  lag  elrcalt  mm  laatallad  to  lagriraa  tka  dynamic  rang  a  of  tka 
3.3  aa A  10.3  m  aavalaagtk  receiver*.  tka  aateaAeA  dynoal c  raaga  amplifiers  at  0.33  aai  1.06  m  par- 
farmed  aa  «ng.cted.  bwovor,  tka  fa  flag  clrealt  far  3.1  sad  10.6  m  mm  dealgaed  for  uaa  with  a  epe- 
elfla  detector.  Thla  detector  fall  ad  Aar  lag  tka  flalA  prograa,  aad  mm  replaced.  Whan  no  ad  atth  tka 
rap  lac ad  dat actor,  tka  gating  latradacad  a  traaalaat  aolaa  that  latarfarad  with  tka  target  return. 
Therefore,  tka  gating  mm  not  need  on  tka  gaaka  Meek  TT/fkOW-TVO  coat a.  gather,  optical  flltere  ware 
need  oa  the  3.1  and  10.6  in  aevelaegth  ayetame  to  prevsat  target  aad  daaaa  aeroeol  eaturatlon  of  re¬ 
ceiver  elaetroalca.  Aa  a  coaaageaaco,  low-level  claar-alr  acattorlng  noma  1  ly  wee  net  okaarrad. 
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(V)  Figure  5  present >  individual  wavelength  beckecetter  signatures  recorded  during  e  light 
snowfall.  A  3d*  optical  filter  was  added  to  the  0.33  \m  end  1.06  ua  wavelength  receivers  to  prevent 
receiver  saturation  by  the  large  amplitude  returns  from  the  enow  near  the  Ildar  van.  Range-resolved 
■now  returns  are  observed  alwoet  to  the  target  with  the  3.8  m  end  10.6  m  Ildar  wavelengths. 


LIGHT  SNOWFALL  IS  JANUARY  1984 


8  bit  -  100  Mz 


100  n/div  3.8  pH  100  H/Div  10.6  p* 


FIOU8C  •  UOAH  tACMCATTM  SKMATURtS  AT  0J1  \M.  3JI  net  tOJ  m  FOB  LIGHT  SNOWFALL 

coNomoNS 


(B)  Figure  6  presents  oscilloscope- photographed  Ildar  beck  scatter  signatures  that  wars  recorded 
near  the  and  of  Trial  310  (2+9  ala).  The  suoke  event  wee  generated  bp  three  white  phosphorous  rockets. 
The  Ildar  returns  free  the  relatively  tenuous  tasks  are  of  greater  amplitude  chon  the  target  return. 
Ales,  the  patchy  snake  structure  results  la  0.33  w  wavelength  tanks  returns  that  am  aaarly  Identical 
la  shape  to  the  target  returns.  Aa  tha  laser  pulse  deration  Increases  (as  for  the  longer  wavelengths), 
the  tasks  patches  ere  spatially  averaged,  aod  the  resulting  make  returns  ere  easily  Identifiable  fron 
target  retene.  (In  addition.  It  eaa  bn  noted  that  tha  gey  elm  recems  la  Flgem  6  are  asm'  difficult 
to  separate  fron  the  snake  returns  for  tha  longer  pels#  eye teas.) 
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TRIAL  310  13  JANUARY  1984  Z  ♦  9  minutes 

8  bit  -  100  mi 


.5? £W: 


wmy 

i  V'/.!  •  ;  . 


100  n/oiv 


.532  jin 


110  m/div 


10.6  jim 


FKJUKES  LtOA*  SACICSCATTIR  SIGNATURES  AT 0J3.  1.0S,3JsnS  IflApai  FOR  WHITE  PHOSPHOROUS 
SMOKE  AT  1*9 MINUTES  OP  T8UM. 310 

(U)  The  four-wavelength  backseat tar  signatures  itorad  on  aagnetlc  tap*  can  be  processed  aa  in¬ 
tensity  updulated  displays  ublch  UXuatrata  tba  tina  and  specs  variability  of  aaoka  and  tar*at  Ildar 
raturaa  occurring  during  a  aaoka  trial.  Figaro  7  praoaata  rack  a  display  for  data  collactad  during 
Trial  218  on  12  January  1984.  Tba  aaoka  ovaet  «aa  ganaratad  by  thraa  rad  pboophoroua  Zunl  rocketu. 
figure  7  lndlcataa  that  tba  aaoka  co^lataly  obacurad  tba  targat  at  all  four  lldar-vsvelragtha .  dovever, 
tba  targat  waa  again  observed  within  a  half-alnuta  at  tba  1.06  wa.  3.8  ua,  and  10.6  ua  wavelengths; 

Tba  targat  was  obacurad  for  about  2  mnutaa  at  0.53  ua.  Detailed  atructura  of  tba  aaoka  dlr'eraioa  la 
observed  at  0.33  and  1.06  \m  but  ant  at  3.8  and  10.6  ua  bacauoa  of  lo^ar  pulaa  langtha  and  lower  de- 
t actor  banduldtha  of  tb>  far  lnfrarad  ayataaa. 


(0)  figure  8  p raa seta  aa  Intensity  nodulated  display  for  Trial  128  conducted  on  16  January  1984. 
free  foreign  pots  was  eery  effective  la  obscuring  tba  targat  at  0.33  and  1.06  ua,  but  not  at 
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CD)  The  first  task  of  ths  analysis  program  was  to  avaluats  tha  response  of  tha  Ildar  racalvar 
from  calibration  data  collactad  during  tha  Smoka  Mask  VX/SBOtf-THO  parlod.  Xecelver  calibration  con¬ 
sisted  of  Inserting  optical  filters  of  known  attenuation  In  front  of  the  detector,  and  of  making  Ildar 
observation* •  Tha  decrease  In  target  returns  Introduced  by  tha  optical  filter  provides  a  measurement  of 
the  response  of  tha  Ildar  racalvar  In  terms  of  tha  number  of  digitiser  counts  par  dB  of  attenuation.  In 
practice,  three  filters  of  different  attenuation  were  used  to  evaluate  linearity  and  response  over  the 
dynamic  range  of  tha  receiver.  The  response  factors  than  are  uaed  to  correct  backscatter  signatures  In 
terms  of  relative  light  input  Into  the  Ildar  receiver.  Pulse- to-pulae  transmit  energy  variations  were 
recorded  and  these  data  were  used  to  normalise  tha  backscatter  signatures  to  a  constant  transmit  level.  ‘ 


TRIAL  m  FOREIGN  SOTS  ts  JANUARY  MM  •  WT/100  MM* 
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FIGURE  •  RANGE/TIME  tNTCNSITY  MOOULATEO  OISPLAY  GENERATED  FROM  FOUR-WAVELENGTH 
LIOAR  BACKSCATTER  SIGNATURES  COLLECTED  DURING  TRIAL  120 

CO)  Lldar  data  collected  during  the  smHBG-l  teats  showed  that  the  target  returns  could  be  In¬ 
terpreted  In  tens  of  transmission  that  agreed  well  with  tranamlaslea*  measured  by  two-ended  trans- 
nteaonetar  systems  (Utha  at  al.,  1983).  For  example,  Figure  9  presents  lldar  and  ML  transmlssooetar 
records  obtained  at  a  wavelength  of  1.08  urn  during  an  IX  screaoer  smoke  event  conducted  on  8  December 
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1982.  Excellent  agreement  1s  obtained,  although  the  two  obsarvational  paths  ware  separated  by  about  3m 
la  the  horizontal  ■  Therefore,  the  Ildar/ car  jet  transmission  values  are  considered  valid,  and  can  be 
used  to  validate  various  aathods  for  Interpreting  aerosol  backacattar  signatures  In  terms  of  optical 
density. 

03)  figure  10  presents  an  example  of  four-wavelength  Ildar- evaluated  transmissions  for  the 
Seeks  Week  VI  Trial  218  (see  Figure  7).  The  snoka  event  was  generated  by  three  rad-phosphorous  rockets. 
The  pre-event  target  returns  are  used  to  normalize  the  transmission  evaluations  to  100Z.  Fluctuations' 
in  the  0.S3  ua  wavelength  any  be  the  result  of  digitiser  errors  for  large  amplitude  (target)  signals. 

Bo waver,  the  fluctuations  are  not  a  problem  at  lower  transmission  values  (e.g.,  smell  target  returns). 
Figure  10  shorn  that  the  attenuation  at  0.33  ua  is  over  two  times  greater  than  at  1.06  ua.  Attenuation 
at  3.8  ua  and  10.6  \m  la  only  about  one-four  of  that  at  1.06  ua.  The  relatively  little  data-polnt 
scatter  in  the  optical  depth  plots  ladlcatos  that  apace  and  time  variations  of  the  obscurant  cloud  did 
not  effect  the  sequential  aultlple-wavelength  laser  transmissions.  Also,  the  wavelength  dependence  re- 
aa Inert  relatively  cosecant  daring  the  test.  Indicating  that  particle  sizes,  shapes,  and  coaposltlona  also 
remain  relatively  constant.  Lower  attenuation  than  expected  at  1.06  ua  for  large  optical  depths  can  be 
explained  by  multiple-scattering,  as  discussed  .later  la  this  paper. 

03)  Analysis  of  Ildar  backscatter  signatures 
In  .terms  of  aerosol  optical  density  along  the 
propagation  path  normally  begins  by  assuming  the 
volume  extinction  coefficient  <r  is  related  to  the 
voluae  backscatter  coefficient  8  by  the  expression 
o  *  kt.  Then,  the  range-corrected  tlngle- 
acacterlng  Ildar  aquation  PC*)  *  Cl(t)  exp 
[-2  / >(t)dX],  can  be  analytically  solved  for  o(X) 
when  k  and  a  boundary  value  of  d  at  range  *o  are 
known  (Klect,  1981;  Vlasee  at  al.,  1969).  Bowever, 

Che  single  scattering  neumptlon  la  not  valid  for 
larga  optical  daptha;  end,  therefore,  tbla  tech¬ 
nique  may  be  of  little  use  for  evaluating  dense 
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FIGURE  9  COMPARISON  OF  TRANSMISSION  HISTORIES 

MCASURCO  8V  TH«  UJSisn  SRI  UOAR  ANO  THg 
NHL  TRAMMtStOMCTgR  FOR  IR  SCREENER  SMOKE 
EVENT  ON  8  DECEMBER  1982 
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OPTICAL  DEPTH  AT  INOICATED  WAVELENGTH  (X)  —  Ux 

FIGURE  10  TRANSMISSION  HISTORY  ANO  OPTICAL  OIPTM  RELATIONSHIPS  FOR  SMOKE  GENERATED 
OURING  TRIAL  21*  (3  RP  ZUNI  ROCKETS! 

ob scar sac*.  The  o  -  kfl  relationship  can  be  Investigated  by  relating  target-evaluated  optical  depth 
0  ■  -laT  •/a dft  to  path-integrated  Ildar  signal /FdK.  For  nail  0,  k  «  a/t  can  be  evaluated  from  a 
calibrated  Ildar.  The  relationship  of.  optical  depth- to- Integrated  Ildar  return  Includes  effects  of  at¬ 
tenuation  and  nultlple'  scattering)  and  therefore,  nay  be  useful  for  Interpreting  single-ended  Ildar 
observations  In  tarns  of  trananlsaloa  of  laser  energy  through  dense  obscurant  events.  However,  previous 
results  (Uthe,  1981)  have  shown  relatively  large  data  point  scatter  between  V  and/Pdlt.  This  nay  be 
a  result  of  experimental  errors  associated  with  Ildar  receivers  of  Halted  dynamic  range  and  bandwidth, 
and  was  a  major  reason  for  extending  the  dynamic  range  and  response  of  the  four-wevelength  Ildar 
receivers  (see  Section  2). 


00  Figure  11  presents  1.06  un  wavelength  .optical  depths  plotted  agsinae  1.06  un  path-integrated 
Ildar  signals  for  data  collected  during  Trial  218.  Although  the  guy  wire  Introduces  additional  uncer¬ 
tainty  in  the  integrated  backscatter' values,  the  data  point  scatter  Is  less  than  that  observed  In 
previous  studies.  This  Indicates  that  Improving  tha  dynamic  range  and  response  of  tha  Ildar  receiver 
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hu  improved  the  relationship  of  cho  Ildar  signal  co 
optical  danalty  tarns.  At  low  optical  dspth  values, 
a  a aar  llnasr  rslatlonshlp  Is  obtained,  indicating 
that  tha  o  •  US  relationship  Is  valid.  At  large 
optical  dapehs,  lass  Ildar  signal  is  observed  be¬ 
cause  of  tha  attenuation. 

(U)  Another  method  Is  being  lnvesclgated  for 
providing  profiles  of  extinction  coefficients.  We 
as  sum  that  the  nsxlnun  target  return  results  frost 
non-scat fared  and  near-forward  scattered  light. 

Other  aultlply-scattered  light  incident  on  the 
target  arrives  later  chan  the  non-scactered  and 
forward-scattered  light.  Therefore,  the  Ildar 
target-derived  extinction  coefficient  (e.g.,  trans¬ 
mission)  Is  representative  of  ms 11-angle  nultlple- 
scatterlng  (i.e.,  the  single-scattering  coefficient 
o(  la  reduced  by  the  forward-scattering): 

»  «  «  -  fa(B)  dm  FIGURE  11  OPTICAL  DEPTM-TO-PATH-INTEGRATED  UOAR  SIGNAL 

J  forward  RELATIONSHIP  OERIVED  FROM  1 .08  pm  WAVELENGTH 

direction  UOAR  SIGNATURES  RECORDED  OURING  TRIAL  218 

If  the  scattering  particles  are  of  snail  sits,  the  forward-scattering  tern  la  snail  and  ammag  for  target 
evaluated  traasalaeion,  even  for  dense  obscurants.  When  tha  scattering  particles  are  large  (*10  ua) ,  the 
forward-scattering  tern  at  shorter  wavelengths  la  large  and  effectively  Increases  lldar-darlved  crans- 
alsslona.  Transnlssoaetsr  valuds  any  .e  iner eased  even  sore  If  receiver  aperture  sites  are  larger  than 
chose  used  for  the  Ildar.  Therefor-,  larger  Ildar-observed  cransnlsslens  are  expected  at  shorter  wave¬ 
lengths  for  obscurants  with  lar-t  particle  slsas.  This  Is  Ulus  traced  in  Figure  12  for  data  collected 
during  a  snowfall  at  the  SAOtf-ORg-B  prograa.  Attenuation  at  10.6  ua  is  about  SOX  greater  (e.g.,  saaller 
trsaaalaalons)  than  attenuation  at  1.06  ua.  Longer  wavelength  systaas  designed  to  better  penetrate  saoka 
obscurants  say  be  less  effective  la  penecraelng  snowfall. 

(U)  The  single-scattering  Ildar  equation  can  be  solved  In  terns  of  the  single  constant  ^  (after 
Ternald  at  al.,  1972): 


where  T  Is  tha  target-evaluated  tranaaiaslon,  end 
the  path- integrated  range-corrected  Ildar  signal  j“* 


It  ' 

Is  the  range  to  the  target.  Then  £  can  be  used  with 
FOOdl,  and  the  solution  aquation,  co  evaluate 
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OPTICAL  DEPTH  AT  INDICATED  WAVELENGTH  (X)  —  Ux 

FIGURE  12  TRANSMISSION  HISTORY  ANO  OPTICAL  DEPTH  RELATIONSHIPS  FOR  LIGHT  SNOWFALL  AT  0100 
ON  11  DECEMBER  1M2 


optical  dtp eh  (or  extinction)  to  rang*  a.  Bacauaa  the  Ildar  cargat  cvaluaead  tranaaisalona  ara  us ad  to 
determine  ,  tha  solution  provides  range- resolved  optical  parameters  that  Include  approximate  aultlpla- 
scattarias  effaces.  An  axaapla  of  rsnge-resolved  optical  depth  and  transmission  at  1.06  urn  Is  presented 
In  figure  13  for  Z+63  seconds  of  Trial  218.  Tha  resulting  extinction  coefficient  profile  la  presented 
In  Figure  14.  The  Ildar  solution  greatly  aapllflad  the  optical  density  of  the  far  side  of  the  obscurant 
cloud,  as  viewed  from  the  Ildar. 


(0)  Another  objective  of  this  study  was  to  investigate  methods  of  describing  spatial  and  temporal  } 

variability  of  obscurant  densities  along  propagation  paths.  One  asthod  Is  Illustrated  in  Figure  13. 

Tha  Ildar  transmission  record  presented  la  Figure  9  (IX  Scraenar  Smoke  —  1.06  ua  wavelength)  was  used 

'  i 

to  coapute  tha  number  and  duration  of  smoke  'holes.'  A  hole  occurs  when  the  transmission  becosMS  greater 
than  a  specified  (threshold)  transmission  value,  and  then  later  becomes  less  than  the  specified  value. 

Figure  13  shows  that  the  number  of  holts  Is  gras test  ae  30Z  transmission.  However,  duration  of  the  holes  / 

Is  greatssi.  for  lower  transmission  values.  Similar  analysis  has  been  conducted  for  transmission 
'patches.'  These  results  provide  a  measure  of  obscurant  affsctlveness  at  DoD  laser  wavelengths.  Ludwig 
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•t  al.  (1984  —  In  this  proceedings)  are  axaninip;. 

Ildar  cvo-dlaanslonal  cross  sactlons  of  events 

to  dsscrlba  ths  spatial  distribution  of  *»:ks  and 
clsar  air  pascals. 


5.  0>)  CONCLUSIONS 

(U)  Lldar  is  /.ha  bast  aatbod  to  raaotsly  ob¬ 
serve  ths  spatial  and  tcaporal  distribution  of  ob¬ 
scurant  clouds.  Although  Ildar  also  providaa  a 
aatbod  to  raaotsly  svaluaea  aarosol  dsnaltlas  and 
partlda  eharactarlstlcs,  tha  tachnlqua  is  limited 
by  receiver  dynaalc  ranga  and  bandwidth,  and  by  un¬ 
cart  aintlaa  in  various  analysis  techniques.  For 
dansa  obscurants,  wulttple-s cat taring  affacts  ara 


FIGURE  13  RANGE -RESOLVEO  OPTICAL  DEPTH  ANO  TRANS¬ 
MISSION  (t.08  pm)  ALONG  THE  OPTICAL  PATH 
AT  Z  ♦  <3  SECONDS  OF  TRIAL  318 


laportant,  sad  analytical'  solutions  to  tha  slngla  scattaring  Ildar  aquation  normally  ara  not  appropriate 
and  ocher  methods  oust  be  used. 


FIGURE  14  RANGE-RESOLVED  VOLUME  EXTINCTION 

COEFFICIENT (t. Miami  ALONG  THE  OPTICAL 
PATH  AT  Z  ♦  83  SECONDS  OF  TRIAL  318 


FIGURE  15  NUMBER  ANO  OURATION  OF  HOLES  AT  A 

WAVELENGTH  OF  IjO  «m  FOR  AN  I R  SCREENER 
SMOKE  EVENT  OF  SNOW-ONE-S 
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(0)  This  paper  presents  two  aathoda  to  provide  useful  Information  on  tLe  density  of  obscurant 
cloud*.  Hcvevar,  om  of  the  aathoda  rail**  on  tha  uaa  of  a  targat  to  evaluate  tha  extinction-to- 
backacaetar  ratio  and  lidar-calibration  conatant.  Tha  othar  aathod  require*  a  targat  to  dariv*  an  inte¬ 
grated  lidar  signal-to-optical  depth  ralatlonahip,  and  aaauaaa  thia  ralationahip  la  eharactariatic  of 
ti*  aaoka  typo  ao  that  It  can  ba  applied  to  aiallar  aaoka  events.  Both  aathoda  can  b*  uaad  to 
evaluate  optical  danalty  along  tha  propagation  path.  Tha  targat  aathoda  would  ba  particularly  uaaful  in 
tha  caaa  of  a  downward-viewing  airborne  Ildar  where  the  ground  aurfaca  provide*  a  target  return,  and  the 
ay a tea  can  mk*  obaarvatlona  over  large  regional  area*,  or  track  aaoka  cloud*  to  greater  downwind  dis- 
tancaa.  Thia  technique  already  haa  been  uaad  to  aak*  vertical  tranaaiaaion  gaaaureaent*  of  foreat  fire 
aaoka  at  0.53  ua  and  1.06  ua  (Uthe  at  al.,  1982).  Of  couraa,  if  jr  i*  daterainad  to  be  a  conatant, 
single-end ad  aaaauraaenta  of  optical  danalty  era  faaaibla. 

(B)  The  data  analysis  prograa  1*  continuing  with  support  froa  the  Aray.  Research  Office,  Geosciences 
Division,  and  final  result*  will  be  presented  in  the  near  future.  Result*  already  obtained  have  shown 
the  laportaace  of  lidar  receiver  dynamic  range  and  response  for  aeking  quantitative  danalty  aaaauraaants 
along  propagation  paths. 
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Appendix  D 

"SNOW-ONE-B  MULTIPLE-WAVELENGTH  LIDAR  MEASUREMENTS' 

Presented  at  SNOW  SYMPOSIUM  III 
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(U)  SNOW-ONE-B  MULTIPLE-WAVELENGTH 
USAS  MEASUREMENTS 


Edvard  E.  Utha,  SRI  International 
Stephen  A.  DaLatcur ,  SRI  International 
John  Livingston.  SRI  International 
Norman  B.  Nielsen,  SRI  International 


(U)  ABSTRACT 

'  (U)  Lldar  observations  of  natural 

and  generated  obscurants  were  made  during 
the  SNOW-ONE-B  field  program.  Laser 
pulses  at  infrared  wavelengths  of  1. Oh 
w,  3.8  ias,  and  10.6  ua  were  transmitted 
along  a  common  path  terminated  by  a  pas¬ 
sive  reflector  located  about  600  m  from 
the  Ildar.  Each  pulse  set  (30  sets/ 
minute)  was  transmitted  within  a  150-us 
time  interval.  Laser  energy  back- 
scattered  by  aerosols  and  reflected  by 
Che  target  was  collected  within  a  12-inch 
telescope  and  optically  directed  to 
appropriate  detectors.  Resulting  range- 
dependent  Ildar  signals  were  logarith¬ 
mically  amplified,  digitised  at  a  3-n 
resolution,  and  stored  on  magnetic  tape 
for  subsequent  analysis. 

(U)  Analysis  examples  are  presented 
illustrating  the  evaluation  of  the  time 
history  of  transmission  based  on  observed 
target  returns.  High  correlations 
between  optical  depths  at  different  wave¬ 
lengths  indicate  that  time  and  space 
variation  of  aerosol  concentrations  heva 
little  effect  on  a  pulse  set.  Linear 
relationships  t>‘-rween  the  optical  depths 
indicate  that  multiple  scattering  effects 
are  minimal.  Slopes  of  the  linear  rela¬ 
tionships  provide  estimates  of  the  rela¬ 
tive  effectiveness  of  the  obscurant  sc 
.each  wavelength. 

(U)  Lldar  transmission  histories  are 
compared  to  transalsoonater  values,  and 
are  analysed  la  terns  of  the  total  tine 
for  which  transmission  was  less  than  any 
specified  value.  The  data  is  described 
la  terms  of  traamlaalon  holes  end 


patches.  In  addition,  preliminary  back- 

scatter  analyses  are  presented. 

(u>  introduction 

(U)  SNOW-ONE- B  was  part  of  a  field 
test  series  conducted  by  the  U.S.  Army 
Cold  Regions  Research  and  Engineering 
Laboratory  (CREEL)  to  determine  the 
effects  of  natural  and  man-made  obscur¬ 
ants  on  ultraviolet  through  millimeter- 
weye  propagation  in  winter  environments. 
The  SNOW-ONE-B  tests  were  conducted  at 
•  Cray ling,  Michigan,  during  November  and 
December  1983.  The  tests  included  elec¬ 
tromagnetic  propagation,  aerosol  charac¬ 
terisation,  and  meteorological  data 
collection. 

(U)  Under  sponsorship  of  the  U.S. 
Army  Research  Office  (ARO),  SRI  Inter¬ 
national  participated  in  the  SNOW-ONE-B 
field  program  by  making  multiple- 
wavelength  Ildar  measurements  of  back- 
scatter  and  attenuation  of  aerosol  clouds 
at  laaer  wavelengths  of  1.06  un,  3.8  un, 
and  10-6  wm.  This  paper  presents  first 
results  from  analysis  of  the  SNOW-ONE-B 
lldar  records.  Turn  analysis  program  is 
continuing. 


(U)  MULTIWAVELENCTH  LIDAR  SYSTEM 

(U)  The  SN0W-0WE-B  experiment  was 
the  first  application  of  a  new  multiple- 
wavelength  van-mounted  lldar  system 
especially  designed  for  remote  character¬ 
isation  of  aerosol  densities  and  particle 
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character  is cic».  «J  shown  in  Figure  1, 
the  lidar  can  transmit  laser  pulses  at 
four  wavelengths:  0.53  urn  and  1.06  urn 
using  a  Nd:7AG  laser  that  Includes  a 
doubling  crystal,  3.8  ua  using  a  DF 
l*»«r,  and  10.6  ua  using  a  CO 2  laser. 

The  lasers  pulse  at  each  wavelength 
sequentially  and  transait  into  the  atmos¬ 
phere  within  a  150-us  interval.  This  Is 


achieved  by  double  pulsing  the  Nd:YAG 
laser  during  a  single  flash- lamp  excita¬ 
tion.  The  laser  pulses  are  directed 
along  a  path  coaxial  with  a  12- in 
Newtonian  telescope.  Collected  back- 
scattered  energy  is  directed  to  an  appro¬ 
priate  detector,  and  the  detector  outputs 
are  logarithmically  amplified.  A  multi¬ 
plexer  circuit  provides  a  single-channel 


waveform  consisting  of  th«  output  of  four 
energy  monitors  (pyroelectric  detector*) 
in  Addition  to  the  four  backscatter  sig- 
natures.  Therefore,  a  single  transient 
recorder  (Biomat  ion  8100)  can  be  used  to 
digitise  the  Ildar  records.  A  clock-gat* 
program  unit  provides  for  transient 
recorder  interrupts.  Data  for  each  pulse 
set  is  contained  in  the  2048  8-bit  Bioma- 
tlon  memory.  A  sample  rat*  of  20  ns  (3-m 
resolution)  is  typically  used,  and  30 
puls*  sets  can  be  transmitted  each  minute 
of  operation.  The  digital  record  for 


each  pulse  set  is  transferred  to  a 
Digital  Equipment  Corp.  LSI-11  micro¬ 
computer,  and  written  onto  magnetic  tape 
(nine-track,  1600  BPI)  along  with  date 
and  time  information.  The  lidar  system 
is  installed  in  a  7-m  van  (Figure  2)  to 
facilitate  the  experimental  setup  at 
remote  field  sites. 

(D)  EXPERIMENTAL  PROGRAM 

(0)  The  SNOW-ONE-B  field  program  was 
conducted  over  flat  terrain  near  Grayling 


Michigan,  during  November  and  December 
1982.  The  lidar  van  was  positioned 
alongside  eransaissometer  Instrumenta¬ 
tion  vans  operated  by  the  U.S.  Navy  and 
the  U.S.  Army  Atmospheric  Science  Labora¬ 
tory.  Ttye  lidar  and  transaissometers 
vere  aligned  to  observe  parallel  paths 
extending  over  a  600-m  range.  A  passive 
(plywood)  target  of  64  square  feet  was 
erected  at  the  far  end  of  the  lidar  path, 
and  transaissomecer  energy  sources  were 
located  near  the  lidar  target.  The  tar¬ 
get  provided  a  constant  reflectivity 
surface  such  that  variation  in  lidar- 
observed  target  returns  can  be  inter¬ 
preted  in  terms  of  atmospheric  attenua¬ 
tion  over  the  lidar— to— target  path. 

(U)  Because  of  eye-safety  require¬ 
ments,  the  0.53  pm  laser  wavelength  was 
not  transmitted  during  the  experiment. 
Therefore,  lidar  observations  were  made 
only  at  the  infrared  wavelengths  of 


1.06  urn,  3.8  um,  and  10.6  urn.  The  lidar 
was  allowed  to  operate  the  first  20  min¬ 
utes  of  each  hour  during  measurement 
periods.  Occasionally,  vertical  lidar 
measurements  were  made  by  placing  a  45° 
mirror  outside  the  laser  van. 

(U)  The  primary  objective  of  the 
field  program  was  to  characterize  elec¬ 
tromagnetic  propagation  through  dense 
snowfall.  Unfortunately,  only  light 
snowfall  occurred  during  the  experimental 
period.  However,  cm  several  days  s  rela¬ 
tively  dense  fog  occurred,  end  measure¬ 
ments  were  made.  (Observations  also  ware 
conducted  for  several  genera  ted- smoke 
events.)  Table  1  presents  a  summary  of 
the  lidar  operations. 

(U)  EXPERIMENTAL  RESULTS  (Preliminary) 

(U)  In  this  section,  we  present  some 
of  the  results  obtained  thus  far  by  s 


Table  1.  Lidar  Maaaurianta  taken  During  SWV-ORE-B  Experiment 


SNOW-ONE-B 

Sate  Collection  Rues 
Grayling ,  Michigan 


Data 

■aa 

■Sul 

mm 

■iiJifl 

Veather 

Conditions 

Data  Type 

E3 

KM 

11/18/82 

17:23 

17:3* 

clear  sir 

csllb 

A 

11/19/82 

16:37 

17:02 

high  humidity 

ealib 

8 

11/29/82 

16:40 

17:02 

high  humidity 

celib 

l 

11/30/82 

7:29 

7:33 

elssr  air 

ealib 

2  , 

12/1/82 

7:22 

7:33 

fog 

horls 

3 

12/1/82 

7:35 

8:01 

f®8 

vert 

3 

12/2/82 

7:14 

7:37 

light  rain 

ealib 

4 

12/2/82 

7:41 

8:02 

light  rain 

Boris 

5 

12/4/82 

16:17 

16:37 

clear  sir 

ealib 

6 

12/4/82 

16:38 

16:43 

clear  air 

ealib 

6 

12/6/82 

9:00 

.9:1* 

elaar  air 

eal&b 

7 

12/8/82 

10:00 

10:30 

light  aaowfall 

csllb 

8 

12/8/82 

11:00 

11:10 

elaar  air 

harts 

9 

12/8/82 

12:03 

12:19 

oil  fog  omoko 

horls 

9 

12/8/82 

13:03 

13:20 

I*  screens  r  smoko 

berls 

10 

12/8/82 

14:03 

14:20 

fly  ash  smoke 

beds 

11 

12/8/82 

13:06 

13:20 

diesel  smoko 

berls 

11 

12/10/82 

22:05 

22:20 

light  saovfall 

horls 

12 

12/10/82 

23:04 

23:20 

snowfall 

berls 

13 

12/10/82 

24:00 

00:20 

light  snowfall 

harts 

13 

12/11/82 

1:00 

1:20 

light  snowfall 

berls 

13 

12/11/82 

2:01 

2:20 

light  snowfall 

berls 

14 

12/11/82 

3:00 

3:10 

light  snowfall 

horis 

14 

12/11/82 

3:11 

3:20 

light  snowfall 

vere 

14 

12/12/82 

8:00 

9:20 

light  snowfall 

berls 

13 

12/12/82 

10:06 

10:20 

light  snowfall 

berls 

15 

12/12/82 

11:00 

11:20 

light  snowfall 

berls 

15 

12/12/82 

12:00 

12:20 

e]iir  sir 

ealib 

15 
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data  analyala  program.  The, analysis  pro¬ 
gram  la  being  directed  to  the  following 
objectives: 

•  Evaluation  of  1.06  ua,  3.8  ua, 
and  10.6  ua  tranaaiaaion  histor- 
iaa  oyar  tha  lidar-to-target  path 
for  aach  obacuration  avant 

a  Evaluation  of  tha  accuracy  of  tha 
lidar  tranaaiaaiona  by  estimating 
thalr  corralation  to  othar  trana¬ 
aiaaion  aaaauraaanta 

•  Derivation  of  tranaaiaaion  dia- 
trlbution  nodal a  that  nay  ba  aora 
uaaful  to  waapona  analyst a  than 
tranaaiaaion  historian 

•  Evaluation  of  tha  relationship  of 
lidar  backacatter  along  tha 
observed  path  to  tranaaiaaion 
aaaauraaanta 

•  Evaluation  of  tha  Information 
content  on  particle  character¬ 
istics  and  aultiple- scat taring 
effects.  In  order  to  form  batter 
estimates  of  absolute  aerosol 
concentrations  from  single-ended 
lidar  aaaauraaanta. 

(U)  Figure  3  presents  a  single 
aultl-vavelength  lidar  backacatter  wave- 
fora  for  an  oil  fog  aerosol  cloud. 

Ke turns  froa  deer  air,  the  aerosol  cloud, 
and  the  solid  target  can  be  identified 
for  aach  infrared  wavelength.  The  solid 


Taanr  Tanner  Tanner 


Figures  Three  Wavelength  Udar Bacfc- 
Scattsr  Signature  for  Oil  Fog 
Fetes—  Along  800m  Feth 
letween  Ltd er  Van  end  Target 

1N2 


target  returns  were  analyzed  for  time 
histories  of  tranaaiaaion  by  assuming  the 
target  is  of  constant  reflectivity  during 
each  experiaent.  Target  returns  during 
dear-air  conditions  were  used  to  estab¬ 
lish  the  100  percent  tranaaiaaion  value, 
and  optical  filters  with  known  attenua¬ 
tion  at  each  lidar  wavelength  were  used 
to  calibrate  the  lidar  receiver  response 
to  target  returns  in  terms  of  atmospheric 
tranaaiaaion. 

(C)  Figure  4  presents  an  example  of 
a  tranaaiaaion  history  result  obtained 
froa  lidar  target  returns  recorded  during 
the  IK  acreener  avent  on  8  December  1982. 
The  event  observed  by  the  lidar  extended 
over  an  8-ainute  Interval,  with  rapid 
changes  in  trausaiaaion  occurring  at  all 
three  infrared  wavelengths.  The  high 
correlation  of  tranaaiaaion  with  time 
between  wavelengths  results  because 
nearly  equal  paths  are  viewed  at  nearly 
equal  times.  Previous  wultiple- 
wavelength  lidar  measurements  Indicated 
that  path  separations— especially  in  the 
vertical— resulted  in  uncorrelated 
transmission  histories  so  that  aultlple- 
wave length  analysis  could  not  be  applied. 

(U)  Also  presented  in  Figure  4  are 
plots  of  the  optical  depth  (U  “  -In  T) 
between  wavelengths.  Again,  the  rela¬ 
tively  good  correlation  indicates 
minimal  time  and  space  variability  of 
transalsslon  between  the  aultiple 
wavelength  observations.  Moreover,  the 
linear  relationships  between  the  optical 
depths  indicate  that  multiple-scattering 
may  not  be  am  important  factor  for  the 
lidar  transmission  analysis.  Tha  plots 
indicate  that  the  IK  acreener  is  17  per¬ 
cent  more  effective  at  10.6  ua  wavelength 
than  at  1.06  ua,  and  20  percent  more 
effective  at  3.8  ua  than  at  1.06  um. 

(0)  Figure  3  presents  a  comparison 
of  lidar  transalsslon  history  at  1.06  ua 
and  the  equivalent  history  recorded  by 
the  3a vy  trsasmlsoometar  system.  Excel¬ 
lent  agreement  is  obtained,  although  the 
two  observat local  paths  were  separated 
by  about  3  to  4  a  in  the  horizontal. 

(0)  Transmission  histories  such  as 
those  presented  in  Figure  4  are  useful 
for  evaluating  tha  effectiveness  of  ad 
aerosol  avent  against  various  electro- 
optical  systems.  However,  models  pre¬ 
senting  transalsslon  values  as  frequency 
distributions  nay  be  more  readily  useful 
to  weapon  system  designers  and  evalua¬ 
tors.  Several  examples  are  presented 
below  to  illustrate  the  type  of  trans¬ 
mission  models  that  can  be  derived  from 
transmission  histories.  Figure  6  pre¬ 
sents  histograms  of  the  total  time 
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vent  on  8  December  1982 
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TRANSMISSION.  T  -  parcant 

"ijur*  8  Histograms  of  tha  Tima  That  Transmission  Was  Lass  Than  a  Spactfiad 
Vaiua,  T,  for  tha  IR  Scraanar  Smoka  Evant  on  8  Oacambar  1982 


lncraaant  that  transmission  was  balow  any 
spaclf lad  vaiua  fpr  tha  ZR  acraanar  saoka 
cast.  For  axaapla,  at  a  wavalangth  of 
10.6  uai,  tha  tr  ansa  is  a  ion  was  lass  than 
50  parcant  for  about  4  ainutaa  during  tha 
saoka  avant.  Tha  transmission  was  lasa 
than  20  parcant;  for  only  about  1  ainuta 
during  tha  saoka  avant.  Thar af or a,  if  a 
waapon  syataa  is  randarad  inoparativa 
only  if  transalas ion  is  lass  than  20  par¬ 
cant  at  tha  10.6  ua  wavalangth,  tha  ZR 
scraanar  saoka  avant  will  ba  affactlva 
for  about  1  ainuta  of  total  tlaa.  How- 
avar,  tha  affactlvc  tlaa  aay  ba  axtsndad 
bacausa  of  vmriabla  transmissions  axttnd- 
iag  ovar  longar  tlaa'  parioda.  Flgura  7 
praaanta  an  analysis,  of  tha  tlaa  duration 
of  transmission  "holaa"  (i.a. ,  tlaa  dura¬ 
tion  whan  transmission  is  graatar  than  or 
aqual  to  a  spaclflad  vaiua).  For 
axaapla,  only  10  holaa  occurrad  for 
eranscission  of  10  parcant,  but  thasa 
can  ba  of  ralativaly  long  duration.  At 
50  parcant,  nor a  holaa  can  Occur  (28  in 
this  axaapla),  but  with  shortar  tlaa 
durations. 

(U)  A  similar  analysis  can  ba  par- 
formad  for  transmission  "patchas"  (i.a., 
tlaa  durations  whan  transmission  is  lass 
than  a  spaclflad  vaiua).  Figura  8  pra- 
santa  a  transmission  "patch"  analysis 
for  tha  10.6  ua  wavalangth  transmission 
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Figura  7  Numbar  and  Tima  Distribution  of  “Holts" 
(Whan  Transmission  is  Graatar  Than  or 
Equal  to  a  Spscifisd  Vaiua,  T)  for  tha  I R 
Scraanar  Smoka  Evant  on  8  Oacambar 
1982 
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Figure  8  Number  and  Time  Duration  of  "Patches" 
(Whan  Transmission  it  Lmc  Than  a 
Spacifiad  Value,  T)  for  the  I R  Seraener 
Smoko  Event  on  8  December  1982 

history  data  presented  in  Figure  4.  For 
patches  with  less  than  50  percent  trans- 
ails  a  Ion,  50  events  occurred  with  a  maxi¬ 
mise  duration  of  30  sec.  Only  11  patches 
with  transnlsslon  less  than  10  percent 
occurred  with  a  maximum  duration  of  only 
6  sec. 

(U)  The  examples  presented  above 
were  derived  fro*  target  returns.  How¬ 
ever,  the  primary  objective  of  the  Ildar 
program  la  to  derive  quantitative  in for¬ 
mat  Ion  using  single-ended  backscstter 
data.  Several  methods  have  been  used  to 
relate  Ildar  atmospheric  backscstter 
observations  to  extinction  or  concentra¬ 
tion.  For  relatively  dense  aerosols, 
single-scattering  formulations  nay  not 
be  valid  and,  therefore,  the  typically- 
used  Ildar  equation  and  its  various  solu¬ 
tion  forms  may  not  be  appropriate. 

(D)  One  method  for  evaluating  atten¬ 
uation  from  backscstter  data  is  to  relate 
the  path-integrated  backscstter  to  the 
optical  depth  evaluated  from  target 
returns.  For  example.  Figure  9  presents 
a  plot  of  these  two  quantities  for  Ildar 
data  collected  at  a  wavelength  of  3.8  urn 
during  the  IX  acreener  test  (c.f.  Figure 
4).  The  Integrated  backscstter  ratio 


was  calculated  using  a  clear  air  return 
obtained  before  the  smoke  release. 

(U)  Although  a  relatively  large 
scatter  of  data  points  occur,  a  near- 
linear  relationship  can  be  evaluated  for 
optical  depths  less  than  1.5  (l.e., 
transmissions  from  100  to  22  percent),. 
This  linear  relation  can  be  used  to  infer 
transmission  from  single-ended  back- 
scatter  measurements.  The  near  constant 
path-integrated  backscstter  and  the  rela¬ 
tively  large  scatter  of  data  points  above 
optical  depths  of  1.5  have  been  assoc¬ 
iated  with  dynamic  range  limitations  of 
the  Ildar  receiver,  and  possibly  with  a 
wide  secondary  pulse  superimposed  o'  .he 
primary  narrow  pulse  of  the  gas  lav*:* 
(this  introduces  false  scatter t  .c 
ranges  greater  than  the  aeror  0.  target) . 

(U)  Figure  10  shows  the  transmission 
of  three  laser  wavelengths  during  light 
snowfall  (approximately  .03  grams/m3) 
starting  at  0100  on  11  December.  During 
this  snowstorm,  transmission  is  generally 
decreased  with  Increasing  wavelength. 
Figure  10  also  contains  the  relationships 
of  the  optical  depths  for  this  period. 

In  sumary,  Oio.6  *  °3.8  >  D1.06- 

(D)  Figure  11  presents  the  relation¬ 
ship  between  optical  depth  and  Integrated 
backscatcer  ratio  for  the  3.8  urn  data  of 
Figure  10.  Because  a  clear  air  profile 
was  net  available  during  this  period, 
the  backscatter  ratio  la  normalized  to 
the  minimum  profile  obtained  for  the 
data  sat.  As  with  the  smoke  data  pre¬ 
sented  In  Figure  9,  there  Is  large 
scatter  in  the  data  points.  Because  of 
equipment  limitations,  identified  by  the  . 
smoke  and  snow  data  at  S NOW-ONE -B ,  back¬ 
seat  tar  analysis  of  relatively  dunse 
aerosol  events  (e.g.,  smoke  tests)  are 
not  being  pursued.  Bather,  extended- 
range  log  amplifiers  and  time-gated 
detectors  are  being  incorporated  to 
expand  the  system  dynamic  rsnge  for  use 
on  future  tests. 

(U)  This  research  was  sponsored  by 
the  Ceosciences  Division  of  the  U.S. 

Army  Research  Office. 
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Transmission  History  and  Optic*)  Depth  Relationships  for  Light  Snowfall  at  0100  on  1 1  December  1982 
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Appendix  E 

"LIDAR  APPLICATIONS  FOR  OBSCURANT  EVALUATIONS" 
Presented  at  SMOKE/OBSCURANTS  SYMPOSIUM  VII 
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LIDAR  APPLICATIONS  FOR  1BSCURANT  EVALUATIONS  (U) 


Edwar-J  E.  Uth« 

SRI  Intarnatlonfl 

Manlo  Park,  California  94025  U.S.A. 

ABSTRACT  (U) 

(U)  Ildar  provldea  a  unique  nethod  for  obaerving  temporal  and  spatial  distributions  of 
obscurants  along  optical  paths.  How aver,  attenuation  and  multiple  scattering  effects  Introduced  by 
very  dense  obscurant  clouds  can  (1)  severely  limit  observational  ranges,  and  (2)  complicate  Interpreta¬ 
tion  of  backscatter  signatures  In  tens  of  absolute  density  and  other  properties.  Keverthelesa,  Ildar 
properly  used  can  provide  needed  Information  that  Is  nearly  Impossible  to  obtain  by  other  means.  This 
paper  reviews  applications  for  obscurant  evaluations  of  several  existing  lidara.  Including  a  four- 
wavelength  van-mounted  and  a  two-wavelength  airborne  system.  Examples  are  presented  that  Illustrate 
the  evaluation  of: 

e  Aerosol  plume  transport  over  complex  terrain  and  water  aurfsces 
e  Cross-plume  diffusion  during  active  convection 

e  Multiwavelength  transmission  j  over  lidar-to-target  horizontal  paths 
e  Vertical -path  transmissions  using  the  earth's  surface  as  a  passive  reflector 
e  Particle-size  distribution  parameters  derived,  from  multiwavelength  data 
e  Backscatter  properties  of  natural  and  senarated  obscurants 
e  Optical  and  physical  density  Inferences  from  single-ended  observations. 

1.  (U)  LIDAR  StSTSfS 

(U)  Lidar  applies  shorr-wavelength  laser  sources  in  radar  fashion  for  remote  observation  of 
scattering  aaroaola  and  reflective  surfaces.  High-speed  digitization  and  logarithmic  protesting, 
coupled  with  computer-bated  video  displays,  provide  for  reel-time  generation  of  pictorial  views  of 
atmospheric  structure.  This  la  derived  from  range-dependent  backscatter  signatures  collected  vltb  high- 
pulse  rate  lasers.  Van-mounted  and  airborne  lidar  systems  have 'been  developed  for  observing  atmospheric 
aerosols  over  large  volumes  with  high  spatial  and  temporal  resolutions.  This  paper  presents  examples 
derived  from  several  lidar  systems  operated  by  SRI.  Although  these  examples  mostly  have  been  derived  on 
environmental  studies,  they  clearly  relate  to  applications  for  obscurant  evaluations.  Syatem  details 
have  been  presented  elsewhere. 

2.  (U)  AEROSOL  PLUME  TRANSPORT  OVER  COMPLEX  TERRAIN  AND  HATER  SURFACES 


(U)  Downward-viewing  airborne  Ildar  systems  era  particularly  wall  suited  for  evaluating  terrain 
effects  on  the  transport  and  diffusion  of  aerosol  plumes.  For  axample.  Figure  1  presents  plume  cross 
sections  collected  with  the  ALPHA-1  (Airborne  Lidar  Plume  and  Haze  Analyser)  at  different  downwind  dis¬ 
tances  from  a  coal- burning  power  plant.  The  plume  was  being  channeled  by  complex  terrain  to  the  west 
of  tha  plant.  Other  data  collected  show  mountain  restrict  Iona  and  barriers  to  plume  transport 
(Uthe,  at  el.,  1980). 

Paper  to  be  Included  In  the  P?  oc -•dings  of  Smoke/Obicurent'  Symposium  VII,  Harry  Diamond  Laboratories 
Adalphl,  Maryland,  26-28  Aprl  1983 
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(U)  Airborne  systens  arc  also  well  suited 
for  aa Icing  observations  over  water  surfaces.  For 
^  exaaple.  Figure  2  presents  an  exaaple  of  air  pol¬ 

lution  layers  over  Los  Angeles,  and  extending 
over  the  ocean  area.  The  data  show  increasing 
density  and  depth  of  the  pollution  layer  during 
the  uiy,  and  indications  of  elevated  wind  flov 
reverse  .ear  the  San  Gabriel  Mountains.  The 

00  suu rr 

'■  t~ -  ALPHA-!  trsnsaits  energy  simultaneously  at  0.53 

n  - 

4Q^  1.06  yk.  The  ratio  of  surface  returns  at 
-m**  these  tmvelangtha  provider  a  good  indication  of 
|m|  land  or  tater  surface,  as  show*  by  Figure  1. 

Tbcnfore,  the  Alfm- 1  can  be  used  to  nap  ter- 
°* <>VA,fr  r,ln  si svat loo  and  surface  type,  to  correlate 

- -  with  observed  atnospherlc  transport  and  diffusion 

- - ,  characteristics. 
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«eaun.uOeCk>«ss  »|«HI  c*m»o*  cotowoo  sivss  (0)  CXOSS-PUME  DITTOS  I  OK  DUS1KC 

ACTIVE  COKVECTIOK 

»*Olus  i  it*  AineouMS  uoan  osssuvanosm  os  tms  rnosssLUMS 

tmucruns  os  a  suevu^Lt  souse  slant  slums  in\  _ l  w-  .  _ 

AT  OisssncNT  oosunuuo  0>st*ncss  snoM  ms  slant  (U)  rs*ur*  4  •nenpies  of  ALPHA-1 

flights  across  a  power  pluue  at  10  ku  downwind 

froa  a  source  located  la  flat  terrain.  Early 

aoralag  observe  does  show  a  well-defined  pluue  above  the  developing  nixing  layer  while  subeegusnt  passes 

scroee  the  pluue  show  structure  of  claar-atr  eonvoctlvo  colls  sad  ontrel - of  the  pluue  particulates. 

Theee  data  show  that  big  he  at  pluue  concent  ret  loee  during  ectlve  convection  are  aeascieted  with  coupenea- 
tlag  douadreft  met  lone.  Darker  araaa  at  the  tap  ef  convective  celle  prehebly  result  by  hygroscopic 
pertuie  growth  at  the  relatively  warn  end  tumid  elr  ie  lifted  use  sealer  ourromedUge.  Theee  date 
illustrate  the  drauecu  sffset  ef  dayttua  ceevect tea  no  aureaal  ceucaatrat loa  dietr ihut lone. 

*•  (0)  wn.TWAmiJCT*  tumshiuums  ovn  lidas-to-taiket  pat** 

<U>  A  fowr-mvelemgth  >Mmu.<H  Ildar  eyeseu  recently  he.  bee.  dovolopod  te  invent igat*  tech- 
nlguaa  far  rauete  aaaaurauaat  ef  unul  apt  teal  sad  phyalcal  daealtlae  and  af  particle  characteristics. 
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LOS  ANGELES.  CALIFORNIA  -  16  DECEMBER  1979 

w  1120 -mo  nr.  i.mw«.  ito*  tnut  cotwe 


Tha  ayataa  crannies  anargy  sc  0.53  -jm  sod  1.06  m  uavslangtha  using  s  naodymlua-tag  Isssr,  3.8  m  using 
a  OF  Isssr  sad  10.6  n  using  s  CO^  Isssr .  Pulsss  from  Cha  thraa  Isaacs  ara  transmitted  coaxially  with  s 
12-inch  Save  on  lan  telescope.  Four-vavalsogth  backseat  ter  and  transmission  ass  suras  ancs  sra  sad  a  along 
cm  saaa  visaing  pach  within  150  us  to  raducs  offsets  of  t lae-varying  sarosol  distributions. 


(U)  A  local  ease  of  ths  four-utvalangth  Ildar  was  conduct  ad  to  avaluata  systan  par  fo  manes  in 
tn~ Irooaancs  stellar  to  thosa  axperlsocad  during  prov.lous  Ildar  operations  at  Aray  tsoga  fscllltlas 
CJtha,  1961).  Ths  Ildar  van  was  locatad 


about  600  a  from  an  8-te  s guars  pasalva 
raflactor  const  rue  t  ad  oa  t  ha  Scaaford  Uni- 
varsity  flald  sits,  both  dust  and  sacks 
aarosols  wars  ganaratad  naar  t  ha  a  id  poise 
of  ths  lidar-to-cargac  path.  Figure  3 
shows  an  oscilloscopa  praaaatatiaa  of  s 
siAgla  Ildar  algaacwrs  showing  ataospheru 
hachscattar  sad  targat  returns  at  each  of 


MTaacs 


c ►a  four  wsvalaagtha.  Ths  largar  ratio  of 
duat-to-«laar-alr  backseat  tar  for  loagsr 


kMcmai  ns  Aano  0*  sunsac*  omasa*  ar  •  sa-ws  awooshus 
smvsunotms  oassavro  av  acaaa  i  Ouaivd  a 
ruOMT  Ovsa  sssTta  ro  saao 
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uevelangchs  is  expected,  as  the  Urge-part  ids  dust  scattering  is  nearly  wavelength  independent.  How¬ 
ever  the  clear -air  scattering  is  proportional  to  the  Inverse  fourth  power  of  wavelength  (assuming  Ray¬ 
leigh  scattering). 


(U)  The  target  returns  provide  a  Seans  to  evaluate  obscurant  tranaalsaion  along  the  lldar-to- 
target  path,  quantitatively,  figure  6  presents  a  tine  history  of  cranseission  derived  from  data  collec¬ 
ted  during  a  rad-enoka  teat.  These  data  indicate  that  the  optical  properties  of  anoks  are  nearly  equal 
at  0.3)  and  1.0*  m,  with  substantially  laaa  attenuation  at  ).l  la  and  alaoat  no  attenuation  at  10.6  us. 
Traa eels aloe  values  evaluated  f roe  target  return*  can  be  used  to  develop  and  validate  txavaaiaslon  values 
derived  froe  the  range-resolved  backseat  tar  data,  so  that  naeauranant*  can  be  nods  along  vertical  and 
slant  paths. 

(U)  The  four-wavelength  Ildar  was  recently  used  to  collact  bocfcacatter  and  traaanlsslon  data  during 
the  SNOU-ont-g  Teats  conducted  at  Grayling,  Michigan;  the  data  are  bow  being  analysed. 

* 
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TARGET  TARGET  TARGE f  TARGET 

043  0m  1  OS  »m  31  104  w* 

I  OUST  /  DU4T  /  DUET  / 

1  IDE  /  ll«  j  104  A*"  / 


104«m  )4«ot 


FIGURE  S  FOUR-WAVELENGTH  HOAR  iACICSCATTER 
SIGNATURE  FOR  DUST  GENERATEO  ALONG 
BOOm  RATH  BETWEEN  LIOAR  AND  TARGET 
Sun  lord  Field  Site.  7  July  IMZ. 


5.  (U)  VERTICAL-PATH  TRANSMISSIONS  USING 

THE  EARTH'S  SURFACE  AS  A  PASSIVE  REFLECTOR 


(U)  During  4  flight  tent  of  the  ALPHA-1,  obser¬ 
vations  were  ud<  downwind  of  a  foreat  flra  located 
near  the  California  coaat  (Uthe  at  al.,  1982).  Visu¬ 
ally,  the  fire  appeared  to  be  contained  within  a 
aall  area;  the  resulting  sacks  pluae  towered  over 
the  source,  but  was  transported  downwind  et  lower 
altitudes. 


(U)  Figure  7  show*  pluae  cross-sections  de¬ 
rived  froa  Infrared  end  visible  baekacattar  signa¬ 
tures  recorded  during  the  second  pass  of  ALPHA-1 
across  tha  pluae.  For  theee  data,  receiver  gain  was 
reduced  so  that  surface  returns  did  not  saturate 


receiver  electronic*;  ae  a  result,  clear  air  bate  layers  wore  not  as  well  observed  on  this  data. 

Cr eater  pluae  attaouatloa  of  the  visible  energy  then  of  tbo  infrared  energy  la  evident  by  tha  absence  of 


pluae  visible  baekacattar  following  penetration  by  the  laser  pulse  of  tha  daneer  pluae  elaaants.  Pluae 
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transmissions  in  the  vertical  ware  evaluated  by  normalizing  to  100  percanc-obaervad  aurface  returna  In 


the  abaence  of  plume  backscacter,  and  by  aeeiaing  conatant  aurface  reflectivity  and  laser- transmit ted 
peak  power  for  each  laser  firing  made  along  the  cross-plume  path.  Lidar  responae  Information  required 
for  quantitative  analysis  of  backscatter  signatures  was  derived  using  standard  calibration  techniques 
with  neutral  density  filters  of  known  attenuation.  Vertical  transmissions  for  one-way  passage  of  the 
laser  energy  through  the  smoke  plume  are  plotted  in  Figure  7.  The  results  show  minimum  plume  trans¬ 
missions  of  about  SO  percent  at  1.06  uo  and 
about  6  percent  at  0.53  us.  The  strong 

.  _  wavelength  dependence  of  attenuation  sug- 


4 

I 

3 


< 


geata  submicron  particle  sixes,  although 
strong  absorption  in  the  visible,  and  weak 
absorption  In  the  Infrared,  by  plume  con¬ 
stituents  could  also  explain  the  observa¬ 
tions.  A  mean  particle  size  Is  derived 
from  the  observations  la  the  following 
section  of  this  paper. 


osmtvto  »«ow  tune  ace  ammns.AT  »as-  eas-am 
mavtLSHSnm  uwe  ns  m. 


6.  (B)  PARTICLE-SIZE  DISTRIBUTION 

PARAMETERS  DERIVED  FROM  MULTIWAVELWGTH  DATA 

(U)  TO  Investigate  che  use  of  multi- 
wavelength  data  for  characterizing  particle 
size  distributions,  an  experiment  was  con¬ 
ducted  using  tranamlssometers  operating  at 
16  different  wavelengths,  and  mounted  across 
a  10-n  long  aerosol  tunnel  facility  located 
at  SRI  International,  as  shown  In  Figure  8 
(Utbe,  1962).  The  tunnel  was  designed  with 
open  ends  for  evaluating  Ildar  techniques. 
Aerosols  consltad  of  particulate  material 
of  different  else  frac irons  and  single-size 
fracltons  of  sin  other  materials.  Particle 
size  evaluations  were  node  by  aultl-stage 
lapse  tor,  sad  by  dir  permeability  saelysls 
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of  packed  powder.  Examples  of  multiple- 
wavelength  extinction  coefficient  data  is  given 
in  Table  1. 

(U)  The  data  base  indicates  that  mean  par¬ 
ticle  size  smaller  then  1-pa  diameter  could  be 
eeti*ac*d.„u*efully  from  aerosol  extinction 
measurements  using  a  single-laser  lidar  system 
operating  at  1.06-  and  0.53-ua  wavelengths 
(Figure  9).  For  larger  mean  particle  sizes  the  noimg*.  (ui  ouumam os Cxscrimcntal paucity 

extinction  ratio  is  near  unity,  and  longer 

wavelength  systems  are  required.  The  data  Indicate  that  a  two- laser  Ildar  operating  at  10.6  and  0.53  urn 
could  provide  estimates  of  mean  particle  size  to  diameters  of  at  least  6  pm. 

(U)  The  two-wavelength  Ildar  data  collected  with  ALPHA-1  (Figure  7)  can  be  used  with  the  labora¬ 
tory  results  presented  in  Figure  9  to  estimate  mean  particle  size  of  the  smoke  plune.  Assuming  the 
wavelength  difference  results  entirely  from  particle  size  effects,  the  laboratory  data  indicate  particle 
diameters  of  about  0.1  pm,  which  is  consistent  with  in-situ  measurements  reported  by  various  investiga¬ 
tors  in  the  literature.  Inference  of  particle  size  is  important  to  evaluate  lidar  backacatter 
quantitatively  in  terms  of  aerosol  extinction  and  concentration. 
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7.  (U)  BACKS CATT EX  PROPERTIES  OF 

NATURAL  AMS  GENERATES  OBSCURANTS 

(U)  As  shown  above,  single-wavelength  Ildar 
systems  can  ba  uaad  to  obaerva  aaroaoi  dlstrlbu- 
tlona  with  high  apaeial  and  temporal  resolution. 
However ,  quantitative  evaluation  of  optical 
(extinction)  or  physical  (concentration)  denai- 
tlea  f roe  backacattar  records  normally  require 
information  or  assumptions  on  the  size,  shape, 
and  composition  of  the  scattering  particles. 
Multiple-wavelength  ays tame  may  provide  suffic¬ 
ient  Information  on  particle  characteristics 
for  evaluation  of  aaroaoi  denaltlea  from  back- 
scatter  data.  The  axaaplea  presented  below 
Illustrate  the  sensitivity  of  Ildar  signatures 
to  particle  characteristics,  and  suggeat  anal  y  a  la 
techniques  for  density  evaluation. 

(U)  Figure  10  preaenta  data  derived  from 
Ildar  observatlona  of  both  black  and  white  smoke 
generated  by  using  the  seme  technique,  but  using 
materials  of  different  chemical  coupoeltlon.  As 
the  density  of  the  plume  was  Increased,  Its  opa-  1 
city  was  measured  by  analysing  the  clear-alr 
Ildar  returns  observed  before  and  after  the  plume 
return;  this  type  of  analysis  yields  a  measure  of 
plume  opacity  that  la  oot  affected  by  uncertain¬ 
ties  la  plume  particle  else,  shape,  and  composi¬ 
tion.  la  Figure  10,  opacity  values  derived  la 
this  meaner  are. plotted  versus  simultaneously 
measured  maximum  plume  backseat ter.  These  data 
show  that  for  visible-wavelength  Ildar  systems, 
approximately  3  49  (factor  of  3)  variation  la 


aidunce  tut  nano  os  ixtinction  at  i  nee  -  me  oau  -pi,  wave  - 

UNCTMS  AS  A  FUNCTION  OF  MAN  ItAUTtAI  DIAMTCR 
(IMFACTOR  MASURSMNTSI  FOR  SUOMICftON  FAATICLt* 

OS  MFFMCNT  cousosmowe 


swung  in  tie  Fumm  rtturns  sumo  aoainst  osaemss  oseivro 
FROMNSAR- AMOFAR.NOg  CHAR  AM  RCTUR**. 
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backacatter  can  occur  for  a  given  plume  opacity,  bacauaa  of  variation!  In  plum*  micro-characteristics. 
Clearly,  evaluation  of  extinction  coefficients  seeming  backacatter-to-extlnctlon  ratios  rust  consider 
particle  composition  (refraction  Index). 


(U)  Figure  11  preaents  Ildar  backacatter  data  related  to  cransmissometer  derived  extinction  data 
for  generated  aerosols  of  the  aaae  composition  (fly  ash),  but  different  size  distributions.  These 
results  show  that  for  a  0.7-ue  wavelength  Ildar  system,  the  backacatter-to-extlnctlon  ratio  is  nearly 


SIOUM  11.  (U)  VOLUW*  BACKSCATTf  a  COCFSICKNTS  (A  ossinvio 
AT  A  UOAA  WAVtLlHOTH  OS  0.7  fan  MLATtO  TO 
WHITS  •  LIGHT  VOLUW*  (XTINCT10N  CO«MICI*NTr 
TDs  iSesHin  IM  WWW  the  bw  IH  »  tlw  ens  at  ttm 
noniimar  aquation  ImiSosWA 


independent  of  the  particle  size.  It  was  also 
shown,  however,  that  at  this  wavelength  the 
backscattar-to-mass  concentration  ratio  is  greatly 
dependent  on  particle  size.  Results  using  a 
1.06~ym  wavelength  Ildar  Indicated  that  the 
backacattar-to-axtlnctlon  ratio  is  more  dependent 
on  particle  size,  but  tha  backseat ter- co-mass  ratio 
is  lass  dependant  on  particle  size  than  at  the 
0.7-um  wavelength,  Later  experiments  have  shown 
even  less  dependence  on  particle  size  for 
backscatter-to-maee  ratios  at  3  to  10  ua  wave¬ 
lengths.  Of  course,  derivation  of  mass  concentra¬ 
tion  from  backacatter  data  requires  knowledge  of 
the  particulate  spclflc  gravity,  and  thus,  volume 
concentrations  nay  be  more  appropriate  for  Ildar 
measurement. 


8.  (0)  OPTICAL  AMD  PHYSICAL  DENSITY  INFERENCES 

FROM  SINGLE- ENDED  OBSERVATIONS 

(0)  The  data  presented  la  Figure  10  illustrated  a  single-ended  method  for  evaluating  the  optical 
density  (opacity)  of  an  aerosol  plume.  This  method  has  been  formally  accepted  by  the  0.S.  Environmental 
Protection  Agency  (Alternate  Method  1  to  Reference  Method  9)  for  objective  measurement  of  visible  emis¬ 
sions  from  stationary  sources.  Experiments  have  been  conducted  to  evaluate  the  application  of  back- 
scatter  to  evaluate  optical  density  (l.a.,  extinction)  as  shown  by  the  data  presented  in  Figure  11. 

Other  data  relating  optical  depth  to  path- Integrated  backacatter  (Figure  12)  shows  that  tone  Informa¬ 
tion  cam  be  obtained,  but  that  dense  aerosol  clouds  can  Introduce  non-linear  relationships  between 
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backscatter  and  extinction;  and  that  vide  acattar  batvaan  thaaa  quantities  can  occur  which  is  probably 
a  result  of  changing  particle  characteristics  within  obscurant  clouds,  nevertheless,  several 
experiments  have  been, conducted  to  evaluate  aerosol  concentrations. 


0  1  2  3  4  8  6  701234  8  6789 

B  (A  -  0.7  pin)  —  km  B  (A  “  10.6  urn)  —  km 
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OVTM  AND  PATH  INTCOMTCO  0ACKSCATT3R  SOS 
A  CLOUD  OS  VIHICUUUt  OUST. 

(U)  To  Investigate  the  use  of  Ildar  for  quantitatively  evaluating  the  emission  rate  cf  fugitive 
particulate  sources,  an  experiment  was  conducted  with  a  controllable  particle  emission  source.  A 
ground-based  Ildar  made  cross-plume  observations  by  scanning  In  elavatlon  at  500  m  downwind  of  the 
smoke  source.  In  addition  to  the  Ildar  observations,  wind  speed  and  direction  were  recorded  at  the 
Ildar  site.  Croes-plOme  backscatter  values  wore  Integrated  and  adjusted  by  wind  speed  and  direction  to 
provide  en  estimate  of  smoke  emission  rats.  The  amcke  generator  was  used  at  three  emission  ratss,  and 
for  each  rate  10-20  pluee  cross  sections  were  recorded.  Figure  13  shows  that  the  resulting  time- 
averaged  croee-pl'tma  backscatter  values  were  linearly  related  to  the  source  amission  rats  for  the 
relatively  low  density  plumes.  These  results  Indicate  the  potential  of  Ildar  for  monitoring  smoke 
plume  particle  concentrations. 

(U)  Analytical  solutions  to  tha  single-scattering  Ildar  equation  have  been  proposed  by  several 
Investigators  as  a  means  to  determine  range-dependent  extinction  coefficients  from  Ildar  signatures  and 
to  correct  for  attenuation  effects  la  concentration  evaluations.  For  example,  Figure  14  presents  a 
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(KlURt  13.  <UI  (ARTICULAT*  (UNION  HATS  EVALUATED  (ROM 

UDAR  (LOTTID  AS  A  (UNCTION  OS  EMISSION  RATS 
8VALUATKO  (ROM  (ARTICLf  (ISDCR. 

Defied  line  i<  the  best-fit  linstr  relation. 


cross-pluaw  sums  concentration  analysis  derived  from  a  ground-based  elevation  scanning  Ildar  systea 
operating  at  a  wavelength  of  .6943  ua-  The  sues  concentrations  were  evaluated  from  the  solution  form: 


where  ■  1/4.34 

*e  *  o/H 

a  «  extinction  coefficient 

M  «  clear  air  concentration 
o 

and  S  ■  range  corrected  logrlthale  Ildar  signature  as  a  function  of  range 

The  paraaeters  Mg  and  were  evaluated  froa  Independently  aaasured  particle  else  distributions  sad  the 
Mia  scattering  theory.  The  result  shows  that  In  this  case,  attenuation  introduced  about  one  order  of 
asgnltude  (10  dS)  decrease  la  the  backseat  ter  fron  the  far  side  of  the  plans.  For  such  danse  pluses, 
•ul tipis  scattering  probably  should  also  be  considered  la  the  analytical  solatlon.  However,  the  croes- 
pli»e  concentration  evaluated  la  the  exaaple  of  Figure  14  agreed  well  with  stack  Mission  data  when 
corrected  for  observed  wind  speed  at  plane  height. 
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9.  (U)  CONCLUSIONS 

(U)  Th«  <xuplu  presented  in  Chi*  paper 
primarily  were  taken  from  environmental  investi¬ 
gation*  of  Ildar  for  quantitative  evaluation  of 
pollutant  distributions  and  densities.  However, 
the  results  illustrate  Ildar  applications  for 
observation  and  evaluation  of  obacurant  clouds. 
Ground-based  scanning  systems  can  nap  cloud 
structure  to  downwind  distances  of  1  to  2  Ins  and 
backscatter  data  can  provide  Information  on 
cloud  optical  and  physical  density.  Longer 
wavelength  (10.6  urn)  systems  are  better  suited 
for  mapping  dense  clouds  because  of  less  attenua¬ 
tion  and  multiple  scattering  effects  (depending 
on  particle  sire  distributions).  Multiple- 
wavelength  Ildar  systems  can  provide  information 
on  particle  characteristics  and  thus  improve 
density  a  valuation*. 


«a 


in 


HOURS  14.  Ill)  EFFECT  os  ATTENUATION  CORRECTION  AND  RESULTANT 
MASS  CONCENTRATION  DISTRIBUTION  (CSNTIRI.  TH* 
CONTOURS  RtSRtSCNT  VALUES  OS  tO  LOQ10  (M/Mo). 
INNSRS  M  IS  THE  SLUMS  MASS  CONCENTRATION  ANO  Vs 
IS  TMS  CLEAR -AIR  RSSSRENCE  MASS  CONCENTRATION 
OS  IOStsm-1. 


<f(V)  v  Airborne  system*  such  as  ALPHA-1  are 
particularly  well  suited  for  obscurant  evalua¬ 
tions  especially  at  distances  greater  then  1  las 
downwind  of  the  source.  The  system  can  map 
terrain  elevations  and  surface  -  type  and  show 

their  effect  on  cloud  transport  and  diffusion.  The  terrain  also  provides  a  reflective  passive  target 
useful  for  evaluating  vertical  transmissions  of  very  dense  aerosol  clouds.  The  two-wavelength 
observations  provide  Information  cn  wean  particle  else  to  help  correct  for  attenuation  and  multiple 
scattering  effects,  and  thus  improve  density  estimates.  A  new  10.6  urn  wavelength  airborne  Ildar 
system' le  being  cousr-ucted  .-at  may  be  the  optima*  apporsen  for  evaluation  of  very  danse  obscurant 
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(U)  The  preparation  of  this  paper  was  supported  by  the  Army  Research  Office  (ARC),  Geosciences 
Division.  The  four-vavelength  Ildar  experiments  also  ara  funded  by  ARQ.  The  ALPHA-1  programs  were 
conducted  for  the  Electric  Power  Research  Institute  and  the  aeroaol  optical  data  were  collected  on 
projects  funded  by  the  U.S.  Environmental  Protection  Agency. 
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